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Abstract. Analogues of stepping-stone models are considered where the site- 
space is continuous, the migration process is a general Markov process, and 
the type-space is infinite. Such processes were defined in previous work of the 
second author by specifying a Feller transition semigroup in terms of expec- 
tations of suitable functionals for systems of coalescing Markov processes. An 
alternative representation is obtained here in terms of a limit of interacting 
particle systems. It is shown that, under a mild condition on the migration 
process, the continuum— sites stepping— stone process has continuous sample 
paths. The case when the migration process is Brownian motion on the circle 
is examined in detail using a duality relation between coalescing and annihilat- 
ing Brownian motion. This duality relation is also used to show that a random 
compact metric space that is naturally associated to an infinite family of coa- 
lescing Brownian motions on the circle has Hausdorff and packing dimension 
both almost surely equal to i and, moreover, this space is capacity equivalent 
to the middle-i Cantor set (and hence also to the Brownian zero set). 



1. Introduction 

Stepping-stone models originally arose in population genetics. The simplest ver- 
sion can be described as follows. There is a a finite or countable collection of sites 
(the site-space) . At each site there is a finite population. Each population is com- 
posed of individuals who can be one of two possible genetic types, say A or B. At 
each site the genetic composition of the population evolves via a continuous-time 
resampling mechanism. Independently of each other, individuals migrate from one 
site to another according to a continuous-time Markov chain (the migration chain) 
on the site-space. 

If the number of individuals at each site becomes large, then, under appropriate 
conditions, the process describing the proportion of individuals of type A at the 
various sites converges to a diffusion limit. This limit can be thought of informally 
as an ensemble of Fisher-Wright dijfusions (one diffusion at each site) that are 
coupled together with a drift determined by the jump rates of the migration chain 



(see, for example, [Shi80|) 



1991 Mathematics Subject Classification. Primary 60K35; Secondary 60G57, 60J60. 

Key words and phrases, coalesce, partition, right process, annihilate, dual, diffusion, exchange- 
able, vector measure, tree, Hausdorff dimension, packing dimension, capacity equivalence, fractal. 

Donnelly supported in part by NSF grant DMS-9505129, and UK EPSRC grants B/AF/1255 
and GR/M14197. 

Evans supported in part by NSF grant DMS-9703845. 

Kurtz supported in part by NSF grant DMS-9626116. 

Zhou supported in part by NSF grant DMS-9703845 and a Michel and Line Loeve Fellowship. 

1 



2 P. DONNELLY, S.N. EVANS, K. FLEISCHMANN, T.G. KURTZ, AND X. ZHOU 



A natural refinement of this two-type diffusion model, considered in [Han90 



DGV95 , is the corresponding infinitely-many-types model. Here the Fisher- Wright 
processes at each site are replaced by mutationless Fleming-Viot processes of evolv- 
ing random probability measures on a suitable uncountable type-space (typically 
the unit interval [0, 1]). 

Much of the research on such interacting Fisher-Wright and Fleming-Viot diffu- 
sion models (see, for example, |BCG86, FG94, FG96, Kle95|) has centred on their 
clustering behaviour in the case when the space of sites is either the integer lattice 
Z'^ or a discrete hierarchical group and the migration chain is a random walk. That 
is, one asks how regions where "most of the populations are mostly of one type" 
grow and interact with each other. The primary tool for analysing this behaviour 
is the duality (in the sense of duality of martingale problems) between these models 



and sytems of delayed coalescing random walks that was first exploited by | ShiS 



One of the factors that lead to interesting clustering is the scaling behaviour of 
the migration process. However, because random walks on in the domain of 
attraction of a stable law and their analogues on discrete hierarchical group only 
have approximate scaling, the role that scaling plays is somewhat obscured. In 
order to make the effect of scaling clearer, related two-type models were considered 
in | EF96| in the hierarchical group setting. In essence, the processes in [ EF96 are 
the result of taking a further limit in which one "stands back" from the site-space 
so that the discrete hierarchical group approaches a continuous one and the random 
walk converges to a "stable" Levy process on the continuous hierarchical group that 
does have exact rescaling. These continuum-sites, two-type stepping-stone models 
have as their state-space the collection of measurable functions x from the site- 
space (that is, the continuous hierarchical group) into [0, 1]. For a state x and site 
e, the value x{e) is interpreted intuitively as the proportion of the population at 
the site that is of type A. 

One of the noteworthy feature of |EF96| is that the limit models are defined by 
specifying moment-like quantities for the associated Feller transition semigroup in 
terms of systems of delayed or instantaneously coalescing Levy processes^ using for- 
mulae that are analogues of the duality relations between the discrete-sites models 
and delayed coalescing random walks mentioned above (see Theorems 3 and 4 of 
|EF96[ ). In particular, the limit models are not defined infinitesimally via a gen- 
erator, SDE/SPDE, or martingale problem formulation analogous to that of the 
discrete-sites models. We note, however, that it should be possible to "stand back" 
in a similar manner from a discrete-sites model where the migrati on chain is sim- 
ple random walk on Z and obtain the process considered in | MT95 ]: this process is 
constructed there as an SPDE on M but is also dual to dela yed coalescing Brownian 
mo tions v ia the same sort of formulae considered in EF96 |. However, the processes 
in [ EF96| that have their semigroups defined in terms of instantaneously coalesc- 
ing Levy processes do not appear to have even a very informal interpretation as 
SPDE-like objects. Rather, a typical value for such a process is a function x such 
that x{e) E {0, 1} for all sites e, and so such p rocess es are more like continuum 
analogues of particle systems (see Theorem 6 of [EF96|). 

The programme of defining continuum-sites models in terms of "duality" for- 
mulae using instantaneously coalescing Markov processes was continued in | Eva97| 
(see Section ^ below for a recapitulation). There the infinitely-many-types case 
was considered and the migration process (that is, the process used to build the 
coalescing system) was taken to be a general Borel right process subject only to 
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a duality condition (with duality here taken in the sense of the general theory of 
Markov processes). Now a state of the process, which we denote from now on by X, 
will be a function x from the site-space into the collection of probability measures 
on an uncountable type-space. For a state x, a site e and a measurable subset G 
of the type-space, the value x{e){G) is interpreted intuitively as the proportion of 
the population at the site possessing types from G. 

We give a more concrete description of the infinitely-many-types, continuum- 
sites process X in Section ^ Under suitable conditions on the migration process, we 
show (at the level of convergence of finite-dimensional distributions) that X is the 
high-density limit of a family oi particle systems with the following description. The 
particles move about in the Cartesian product of the site-space and the type-space. 
The site-space-valued components of the particles evolve according to independent 
copies of the migration process. When particles collide in the site-space a type 
is chosen at random from the types of the particles participating in the collision 
and the types of all the participating particles are changed to this randomly chosen 

type. 

One of the open problems left by |Eva97 | was to determine conditions on the mi- 
gration under which the process X (which is again a Feller process) has continuous 
rather than just cadlag sample paths. In Theorem 7.2, we establish the sufficiency 
of a mild condition to the effect that the coalescing system doesn't coalesce too 
rapidly. The condition holds, for example, for all Levy processes on R. 

By the same argument as in the proof of Proposition 5.1 of |Eva97|, it is pos- 
sible to show that if the migration process is a stable process on the circle T that 
hits points, then, for fixed < > 0, there almost surely exists a random countable 
subset {fci, A:2, . . . } of the type-space such that for Lebesgue almost all e S T the 
probability measure Xtie) is a point mass at one of the ki. That is, rather loosely 
speaking, at each site all individuals in the population have the same type and the 
total number of types seen across all sites is countable. We improve this result in 
Theorem 10.2 for the case of Brownian motion migration on T by showing that 
the total number of types is, in fact, almost surely finite and such a result holds 
simultaneously at all positive times rather than just for fixed times. 

The primary tool used in the proof of Theorem 10.2 is a duality relation be- 
tween systems of coalescing and annihilating Brownian motions that is developed 
in Section ^. This relation enables us to perform detailed computations with the 
system of coalescing Brownian motions that begins with countably many particles 
independently and uniformly distributed on T. 

The latter process is an interesting object in its own right. In particular, it can 
be used to define a random metric on the positive integers by declaring that the 
distance between i and j is the time until the descendents of the i^^ and j*'' particles 
at time zero coalesce. In Theorem 11.2 we adapt the methods of [Eya98| to show 
that the completion of the integers in this metric is almost surely compact, with 
H aus dor ff and packing dimensions both equal to ^. Moreover, this space is capacity 
equivalent to the middle-^ Cantor set (and hence also to the Brownian zero set). 



Notation 1.1. Write N:={1,2,...}. We will adopt the convention throughout that 
the infimum of a subset of M or N is defined to be oo when the subset is empty. 
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2. Coalescing Markov processes and labelled partitions 

Suppose that i? is a Lusin space and that (Z, ) is a Borel right process on E 
with semigroup {Pt}t>o satisfying Ptl = 1, i > 0, so that Z has infinite Hfetime (see 
|Sha88 for a discussion of Lusin spaces and Borel right processes) . Suppose that 



there is another Borel right process {Z,P^) with semigroup {Pt\t>Q and a diffuse, 
Radon measure to 7^ on (S, £) such that for all non-negative Borel functions 
on f,g on E we have J m(de) Ptf{e)g{e) = J m{de) f{e)Ptg{e) (our definition of 
Radon measure is that given in Section III.46 of [ pM7S( |). The space E is the site- 



space and Z is the migration process for the continuum-sites stepping-stone model 
X, whereas Z will serve as the basic motion in the coalescing systems "dual" to X. 

We remark that our assumption on the Markov processes Z and Z is not quite 
the usual notion of weak duality with respect to to (see, for example. Section 9 of 



GS84|); in order for weak duality to hold we would also require that P'"-a.s. (resp. 
P™-a.s.) the left-limit Z{t-) (resp. Z{t-)) exists for aU t > 0. 

The following notation will be convenient for us. Given e — (ei, . . . , e„) S i?", 
for some n e N, with e, ^ for i ^ j, let Z"" = [Z"^ , . . . ,Z'=") be an E"-- 
valued process defined on some probability space (fi, J^, P) such that Z'^^ has the 
distribution of Z under P'^' and Z'^'^ , ■ • • , are independent. 

We now define the system of coalescing Markov processes Tf' associated with 
Tf^ . Adjoin a point f, the cemetery^ to E to form E^^ := E U {f}. Construct the 
(i?''^) "-valued process Z'^ = {Zf, . . . , Z^) inductively as follows. Suppose that times 
To < . . . < Tfc < 00 and sets {1, ... ,n} —: Oq ^ ■ ■ ■ ^ Qk ^ {1} have already 
been defined and that Z'^ has already been defined on [0, Tk[. If Tk = 00, then just 
set Tfc_|_i :— 00 and Qk+i ©fc- Otherwise, put 

(2.1) Tk+i := inf{i > Tk : 3t,j G Qk, i^j. Z'^{t) = Z-'^it)}, 



(2.2) 



©fc, if Tfc+i = 00, 

{i e Gfe < j e Gfc, Z^^iTk+i) = Z^-^{Tk+i)} , otherwise. 



and 

(2.3) Z<r{t) 



Z^^{t),Tk<t<Tk+i, ifieOfe, 
t: Tk <t < Tfc+i, otherwise. 



In other words, the coordinate processes of the coalescing Markov process Z'^ evolve 
as independent copies of Z until they collide. When two or more coordinate pro- 
cesses collide (which happens at one of the times with < < 00), the one with 
the smallest index "lives on" while the other coordinates involved in the collision 
are sent to the cemetery f . The set Qk is the set of coordinates that are still alive 
at time Tk- As the following lemma shows, for TO®"-a.e. e almost surely only one 
coordinate process of Z'^ is sent to the cemetery at a time in the construction of Z'^ . 
(Recall that {Z, P^) is said to be a Hunt process if Z has cadlag sample paths and 
is also quasi-left-continuous; that is, if whenever Ti < T2 < . . . are stopping times 
for Z and T = sup„ r„, then P^{lim„ Z(T„) = Z{T), T < 00} ^ P^{T < 00} for 
ah z e E.) 



CONTINUUM-SITES STEPPING-STONE MODELS 



5 



Lemma 2.1. Let Y be an E-valued Markov process on some probability space 

(f2,jr, P) with the same law as Z under P'^ := q{dz) , where q is a prob- 
ability measure on {E,£) that is absolutely continuous with respect to m. Let 
{T,V) be a [0,oo[x E-valued random variable that is independent of Y. Then 
P{y(r) = V} = 0. Moreover, if Z is Hunt process, then f{Y{T-) = V} = 0, 
also. A similar result holds with Y replaced by a process Y with the same law as Z 
under P" := Jj^q{dz)P^. 

Proof. For fixed t > and v G E we have, writing h for the Radon-Nikodym 
derivative of q with respect to m, 

(2.4) F{Y{t)=v}= [ m{dz)h{z)Ptl[^y{z)= [ m{dz)l[^y{z)Pth{z) = 

Je Je 

by the duahty assumption and the assumption that m is diffuse. Moreover, under 
the Hunt assumption, 

(2.5) F{Y{t) ^ Y{t-)} = 0. 

The result now follows by Fubini. □ 

It will be convenient to embellish Z*^ somewhat and consider an enriched process 

defined below that keeps track of which particles have collided with each other. 

Let n„ denote the set of partitions of N„ := {1, . . . ,n}. That is, an element 
TT of n„ is a collection tt = {Ai, . . . , Ai^} of subsets of N„ with the property that 
\J- Ai = Nn and Ai C] Aj =0 for i ^ j. The sets Ai,. . .Ah are the blocks of the 
partition tt. Equivalently, we can think of n„ as the set of equivalence relations on 
N„ and write i ~7r i if i and j belong to the same block of tt G n„. 

An E -labelled partition of N„ is a collection 

(2.6) X = {{Ai,eAj,... ,{Ah,eA^)}, 

with {Ai,... ,Ah} e n„, {bai,.-. ,eAh} C E, and e^* ^ ba^ for i 7^ j. Let 
A„ denote the set of E'-labelled partitions of N„. Put a(A) := {^1, . . . ,Ah} and 

e{X) := {eA)Aea(x)- 

For e e with a ^ ej for i ^ j,we wish to define a A„-valued process (the 
process of coalescing Markov labelled partitions) with the following intuitive descrip- 
tion. The initial value of is the labelled partition {({!}, ei), . . . , ({n}, e„)}. As 
t increases, the corresponding partition a{(^^{t)) remains unchanged and the labels 
e{('^{t)) evolve as a vector of independent copies of Z until immediately before two 
(or more) such labels coincide. At the time of such a collision, the blocks of the 
partition corresponding to the coincident labels are merged into one block (that is, 
they coalesce). This new block is labelled with the common element of E. The 
evolution then continues in the same way. 

More formally, we will take to be defined in terms of as follows (using 
the ingredients and 8^ that went into the definition of Z^). The corresponding 
partition-valued process := a{('^) is constant on intervals of the form [Tfc,Tfe-)-i[ 
and ^^(ro) := {{!},... ,{n}}. Suppose for fc > that ^^(tq),... ,C^(Tfc) have 
been defined and r^+i < 00. Let ^''(r^.+i) be the partition that is obtained by 
merging for each i e 6fe_|_i those blocks of ^''(r^) whose least elements j are such 
that Z^'{Tk+i) = Z'^'{Tk+i). Thus each block A of ^^(r^+i) is such that the least 
clement min^ of A is the unique element i £ A for which Zfirk+i) ^ f. The 
definition of is completed by labelling each block A of the partition S,'^{t) with 
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For 1 < i < n, put = (7^ , . . . , 7^), where 

(2.7) 7^(0 :=min{j 0, 
and write 

(2.8) T^{t) {7f (0 ■.l<t<n}^{j: ZJ{t) ^ f} 
for the set of surviving indices at time t. Note that T'^{Th) — 8^. 

3. The state-space S of the stepping-stone process 
We need some elementary ideas from the theory of vector measures. A good 



reference is [DU77|. Recah the measure space {E,£,m) introduced in Section 



and let B be a Banach space with norm || • ||. We say that a function (p : E ^ B 
is simple if = X^iLi^^ilfii for xi, . . . , Sfc £ B and Ei,... ,Ek G £ for some 
fc e N. We say that a function (p : E ^ B \a m-measurable if there exists a 
sequence {(f>n}neN of simple functions such that lim„_^oo |l'/>n(e) — 4'{^)\\ — for 
m-a.e. e € E. 

Write K for the compact, metrisable coin-tossing space {0, 1}^ equipped with the 
product topology, and let /C denote the corresponding Borel cr-field. Equivalently, 
/C is the cr-field generated by the cylinder sets. 

Write M{K) for the Banach space of finite signed measures on {K,IC) equipped 
with the total variation norm || • ||M(i<')- Let L°°{m,M{K)) denote the space 
of (equivalence classes of) m-measurable maps ji : E ^ M{K) such that 
esssup{||/x(e)||M(i<-) : e e £'} < 00, and equip L°°{m,M{K)) with the obvious 
norm to make it a Banach space. 

Write C{K) for the Banach space of continuous functions on K equipped 
with the usual supremum norm || • ||c(a:)- Let L^{'m,C{K)), denote the Banach 
space of (equivalence classes of) m-measurable maps fj. : E ^ C{K) such that 
J m^de) \\fJ'{e)\\c{K) < 00, and equip L^{m,C{K)) with the obvious norm to make 
it a Banach space. 



From the discussion at the beginning of §IV.l in [DU77| and the fact that 
M{K) is isometric to the dual space of C{K) under the pairing {v^y) i— > {v^y) — 
/ u{dk)y{k), V e M(X), y e C{K), we see that L°°{m,M{K)) is isometric 
to a closed subspace of the dual of L^{m,C{K)) under the pairing {fi,x) 1-^ 
Jm{de){ii{e),x{e)), /x e L°°{m,M{K)), x e L^{m,C{K)). Write MiiK) for 
the closed subset of M{K) consisting of probability measures, and let S denote 
the closed subset of L°°{m, M{K)) consisting of (equivalence cl asses o f) maps with 



values in Mi{K). From Corollary V.4.3 and Theorem V.5.1 of [ DS58 | we see that, 
as L^{m, C{K)) is separable, S equipped with the relative weak* topology is a com- 
pact, metrisable space. From now on, we always take S to be equipped with the 
relative weak* topology. 

We think of the set K as the space of possible types in the infinitely-many- 
types, continuum-sites, stepping-stone model X we will define in Section |4[ As we 
remarked in Section the type-space for infinitely-many-types models is usually 
taken to be [0,1]. However, from a modelling perspective any uncountable set is 



equally suitable, and, as pointed out in | Eva97| , the set K is technically easier to 



work with. The set E is the corresponding space of sites. The intuitive interpreta- 
tion is that /X g S describes an ensemble of populations at the various sites: fj,{e){L) 
is the "proportion of the population at site e € E that has a type belonging to the 
set L e /C". 
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Remark 3.1. One can think of S as a subset of the space of Radon measures on 
ExK hy identifying ^ £E with the measure that assigns mass m{de) ii{e){B) to 
the set Ax B, where A £ £ and B £ JC. The topology we are using on S is not the 
same as the trace of the usual topology of vague convergence of Radon measures. 
However, the corresponding Borel cr-fields do coincide. In particular, we can think 
of S-valued random variables as random Radon measures on E x K. 

For n e N let M{K") (respectively, C{K^)) denote the Banach space of finite 
signed measures (respectively, continuous functions) on the Cartesian product 
with the usual norm || • \\]\i(K") (respectively, || • ||c(i<'"))- With a slight abuse of 
notation, write ( • , • ) for the pairing between these two spaces. 



Definition 3.2. Given £ ^^(m®", C(K")), define /„( ■ 
of continuous real-valued functions on S) by 



£ C(^) (:= the space 



(3.1) 



™«"(de) ((X) Me.), 0(e)) 



i=l 



) H{e,)idki)(t>{e)ik.), ^ G S. 



i=l 



Write / for h. 



4. Definition of the stepping-stone process X 



Theorem 4.1 below is Theorem 4.1 of |Eva97|. As discussed in Section |l|, it is 
motivated by the characterisation of infinitely-many-types, discrete-sites stepping- 
stone processes via duality with systems of delayed coalescing continuous-time 
Markov chains (see [ pGV95| ] and |Han90( |). Recall that (f^, JP,P) is the probability 
space on which the processes Z'^, Z'^, ■f'^, et cetera are defined. 

Theorem 4.1. There exists a unique, Feller, Markov semigroup {Qt}t>o on S 
such that for all t > 0, fi £ E, (jj £ L^{m®", C{K")), n£N, we have 



Qt{ti,dv)In{v\4') 



(4.1) 



'((ie)P 



^i{z'r{t)){dk,)m{k^m^... , 



jer-(t) 

Consequently, there is a Hunt process, (X, Q^), with state-space 5 and transition 
semigroup {Qt}t>o- 



Remark 4.2. The integrand P[- • •] in (4.1) should be interpreted as on the m*^"- 
nuU set of e such that e,; = Cj for some pair («, j). The integral inside the [ ] is 
over a Cartesian product of copies of K, with the copies indexed by the elements 
of V'it). 



Remark 4.3. The following equivalent formulation of Theorem 4.1 will be useful. 
For n G N let Z^"! = (^i"', . . . , ^IT^) be an iJ"-valued process defined on a a-finite 
measure space 

(rjN^ jc-N^pN)^ with 



(4.2) 



£Ay.^j. 



'{de)¥{Z'' £ A). 



8 P. DONNELLY, S.N. EVANS, K. FLEISCHMANN, T.G. KURTZ, AND X. ZHOU 



Define Zl"!, ^["1, 7["1 and rl"! from in ttie same manner ttiat Z^, 7^= and 



V were defined from Z®. Tlie riglit-liand side of (4.1) is just 



(4.3) 



M(Z"'(t))(rf%)'/'(Z["l(0))(A: 



jerH(t) 



Remark 4.4. As we noted in Remark we can tfiink of the process X as taking 
values in the space of Radon measure on £^ x ii' by identifying Xt with the random 
measure that assigns mass Jj^m{de)Xt{e){B) to the set Ax B, where A G £ and 
S e /C. A standard monotone class argument shows that if 4^ is any non-negative 
Borel function on E^^ x K^, then 



(de) 



|Xt(e,)(dfc,)i/.(e,k) 



(4.4) 



'(rfe)P 



n{Zj{t)){dkj)ip{e, k^f(t), ■■■ , k^iit)) 



^^{z\"\t)){dk,)^P{Z^-^0),k 



jerH(t) 



7r(t)' 



'^l">(t)) 



5. A PARTICLE CONSTRUCTION FOR THE STEPPING-STONE MODEL X 

In this section we first construct a finite particle model in which particles move 
through E X K , where we recall that E is our site-space and K is our type-space. 
The i?-valued components of the particles move independently according to the 
dynamics of the migration process Z. The particles interact only when they are 
located at the same site in E, and the interaction that occurs is that the type of 
one of the particles is replaced by the type of the other. The particle whose type 
"wins" is chosen at random from the two particles, with both outcomes equally 
likely. For our purposes here, we assume that the types are constant except for 
these replacement interactions, although we could allow "mutation" of the types 
between the replacement interactions. 

Under suitable conditions on the migration process, we then pass to a high- 
density limit and obtain a process taking values in the space of Radon measures p 
on E X K with the property that p(A x K) = m{A) for A £ £. Recalling Remark 
3.1, we can think of the limit model as a S-valued process, and we establish that 
as such it has the same finite-dimensional distributions as the continuum-sites 
stepping-stone process X. 

Throughout this section we will work on a probability space (r2,JF, P) and we 
will assume the following hypothesis (the definition of a Hunt process is recalled in 
Section 

Assumption 5.1. The processes Z and Z are Hunt. 

For completeness, we recall the following definition and some of its consequences. 

Definition 5.2. Let {S,S) be a measurable space, and let be a cr-finite measure 
on S. Say that a map N from Ct into the collection of measures on {S, S) is a 
Poisson random measure with mean measure v if 
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a) For each A G S, N{A) is a {0, 1, . . . , oo}-valued random variable. 

b) For each A E S with iy(A) < oo, the random variable N{A) is Poisson dis- 
tributed with parameter i^{A). 

c) For Ai,A2, . . . e «S disjoint, the random variables N{Ai), N{A2), . . . are in- 
dependent. 



Remark 5.3. Assume that u is diffuse. Then for x G S, N{{x}) must be zero or 
one, and so we can identify N with its support. We will write a; e A'' if N{{x}) = 1. 

Note that 



(5.1) 



/ N{dx)f{x) 



.xeJV 



= / iyidx)fix), f e L\iy), 
Je 



and more generally, for / e L^(i/®"), 



(5.2) 



E 



Xl, 
Xi 



... ,Xn&N 

7^ Xj, 



[ l^^"{dK)f{x,,... ,Xn). 

Je^ 



5.1. Finite particle systems. Fix a non-zero diffuse finite measure vq on E 

and a probability kernel fi : E x JC ^ [Oil]- Write D^[0, oo[ for the Skorohod 
space of cadlag i^-valued paths and let M denote a Poisson random measure on 
De[0,oo[xK with mean measure 



(5.3) 



(recall that is the law of Z starting at z € E). Thus the push-forward of M by 
the map {(, k) ^ ^(0) (:= the value of the path Q at time 0) is a Poisson random 
measure on E with mean measure vq. More generally, the push-forward of M by 
the map (C, ^) > C(^) is a Poisson random measure on E with mean measure vt, 
where vtiH) = /g vo{dz) Pt{z, H). We assume that vt is diffuse for each < > 0. By 
our duality assumption, this will certainly be the case if z/q is absolutely continuous 
with respect to m. 

Enumerate the atoms of M as {Zi,k\), . . . ,{Zj,k°j) in such a way that the 
conditional distribution of this collection given J = j is that of j i.i.d. De[0, oo[xK- 
valued random variables with common distribution 



(5.4) 



FxG^ uoiEy^ j Mdz) P'{F) ii{z, G). 



We wish to define a collection ki, . . . , k,/ of AT -valued processes in such way 
that the collection [Zi, ki), . . . , {Zj, kj) has the dynamics described above: that 
is, we think of Ki{t) as the type of the particle Zi at time t, and, after two or 
more such particles collide, the particles participating in the collision must be of 
the same type with the common type selected at random from among the types of 
the participating particles (with each possible outcome equally likely). 

Suppose that on the probability space (fi,^, P) we also have defined for each 
k G N a collection {6ik, i G N} of i.i.d. random variables uniformly distributed on 
[0, 1]. We will implement a specific construction of the 9ik below. Define ki, . . . ,k,j 
and times fo < n < . . . (with f/j < f/j+i when f/j < oo) as follows. 
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Put Ki{0) = and fo = 0. Suppose that fo, . . . have aheady been defined 

and, for 1 < i < J, the processes Ki has been defined on [0, ffe] (or [0, oo[ if ffc = oo). 
If ffc = 00, then the definition of Kj, 1 < i < J, is complete and just define = oo 
for i> k. Suppose, then, that ffc < oo. Put 

(5.5) ffe+i := M{t > ffc : Zi{t) = Zj{t), Ki{Tk) ^ Kjih), some i ^ j}. 

Put Kiit) := Ki{Tk) for Tk < t < ffe+i and 1 < i < J. If ffe+i ~ oo, then this 
completes the definition oi Hi, I < i < J. Otherwise, if ffe+i < oo, then define 
Kj(ffe+i), 1 < j < J, as follows. Let fi(ffe+i) := {j : Zj{fk+i) = -^i(ffc+i)}, 
and let 7i(ffc+i) G fi(ffc+i) satisfy %.(ffe+i),/c+i < ^'j,fe+i for all j € Ti{Tk+i). We 
set K;i(ffe+i) = K^,(fj^^j)(ffc). Note that if Zi{fk+i) = Zj{fk+i), then ri(ffc+i) = 
fj(ffc+i), 7i(ffc+i) = 7j(ffc+i), and K,{fk+i) = Kj{Tk+i). 

Our requirement that the types of colliding particles be changed to a type inde- 
pendently and uniformly selected from those of the participants in a collision will 
be met if for k G N, the collection {9ik, i e N} is independent of V ^f^_^, 
where {^f}t>o is the filtration generated by {Zi, . . . Zj) and {!F^}t>o is the filtra- 
tion generated by {k\, . . . , kj). In particular, the distribution of the process {Z, k) 
will be the same regardless of how we define the {9ik} as long as for each k, the 
conditional distribution of {Oik} given is i.i.d. uniform on [0, 1]. 

Wc note that P-a.s. there exists ^ e N such that ti — oo, so that the above 
construction does indeed lead to a value oi Ki{t), 1 < i < J, for all t >0. To see this, 
let Rh{t) — {I < j < J ■ Kj{t) = Kh{t)} for 1 < /i < J and Q <t < sup^. ffc. Since 
there are only finitely many particles, ]P{i?/t(ffc) C Rh[Tk+i) C • ■ • | Tf^ V J^^^} > 
2~'^ > 0. Consequently, cither there exists f ^ < oo such that (ft) = {1, . . . , J} 
or there exists a time after which Zh does not collide with any particle having a 
diff'erent type. 

Now wc will give an explicit construction of the {9ik} which leads to a useful 
construction of our particle system {Zi,Ki), . . . , {Zj, kj). We assign to each parti- 
cle a distinct [0, l]-valued initial level f/°, 1 < i < J, at time and use these initial 
levels to define a family of [0, 1] valued processes of levels {/7i(t)}t>o, 1 < i < ^• 
The {0ik} will be defined using these level processes. We will assume that the con- 
ditional distribution of given M is that of J i.i.d. random variables uniformly 
distributed on [0, 1]. This assumption implies 

J 

i=l 

is a Poisson random measure with mean measure 

(5.7) FxGxH^ [ uo{dz)P''{F)n{z,G)l{H), 

JE 

where I denotes Lebesgue measure on [0, 1]. We define On Uf . For < t < fi, 
set [/,(<) := C/,(0) := Uf . If fi < oo and |fi(fi)| = 1, then put [/^(fi) := t/i(fi-). 
If fi < 00 and r;(fi) = {zi,... n > 1, then put C/i,(fi) := Ui^^{Ti—) where 

(Tl , . . . , (T„ is a uniform random permutation of 1 , . . . , n selected independently of all 
other quantities. Observe that C/i(fi),... ,C/j(fi) are conditionally i.i.d. uniform 
on [0,1] given Tf^y T'i^. Define 0^ := [/^(fi). Put Ui{t) := C/i(fi), fi < t < f2. 
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We continue inductively, at each time < oo randomly permuting the levels with 
indices in each Ti{Tk) and defining 0i.k+i — Ui{fk). 

Although the level assigned to a particle may change at the time of a collision, 
since these changes only involve the permuation of the assignment of the levels, 
the set of levels is the fixed random setU := {U^}. Consequently, we could index 
the particles and their types by their corresponding level; that is, for u G U, define 
Zu{t) = Zi{t) and Ku{t) = Ki{t) if and only if Ui{t) = u. Since the particle assigned 
to level u changes only when the newly assigned particle is at the same location 
as the previously assigned particle, the strong Markov property implies that the 
processes {Z^, u e hi} are conditionally independent given hi and {Zu{Q), u € h(}^ 
and conditionally each Z^ is a Markov process with transition semigroup {Pt}- 
Note that 

(5.8) ffc+i = M{t > Tk : Zu(t) Zy{t), Ku{t-) ^ Hv{t-), some u ^ v}, 

and if we define f„(ffc) -.^ {v & hl : Zy{fk) = Zuih)} and %{Tk) ■= min(r„(ffe)), 
then Ku(ffc) = Kj^(fk)i^k—)- That is, if two or more particles collide, the particles 
involved in the collision "look down" to the lowest level particle at the same location, 
and change types to the type of that particle. (We note in passing that this construc- 
tion is reminiscent of the "look down" construction of the Moran model in [ DK96| .) 



Consequently, if we start with a Poisson random measure on De[0,oo[xK x [0, 1] 

with mean measure specified by ( ^.T) ), then a particle model 

ueu 

is completely determined by the requirement that whenever two or more particles 
"collide" the types of the higher level particles involved in the collision switches 
to the type of the lowest level particle in the collision. This observation allows 
us to extend the construction to systems with infinitely many particles with mild 
additional assumptions. 

5.2. Particle systems with stationary location processes. We now want to 
extend the construction of the previous section to arrive at a model in which the 
distribution of locations of particles is stationary in time, and we want to allow 
for the poss ibili ty of there being infinitely many particles. We emphasize that 



Assumption 5.1 is in force throughout this section. 



Consider a Poisson random measure 

(5.11) ys 



E"(z„ 



on De[0,oo[xK X [0,1] with mean measure specified by (5.7) with vq — m. By 
the assumption that m is Radon, there exist open sets Ei Q E2 Q ■ ■ ■ such that 
m{En) < 00 for all n € N and E ~ [J^ En (of course, if m{E) is finite we can take 
En^ E for all n G N). Put h(„^ {uGU : Zu{0) G E^}, and note that 



12 P. DONNELLY, S.N. EVANS, K. FLEISCHMANN, T.G. KURTZ, AND X. ZHOU 



is a Poisson random measure on De[0,oo[xK x [0,1] with finite mean measure 
specified by (5.7) with — m(- n E„). As in the previous subsection, we can 
construct a corresponding finite particle model 

(5-13) *N(i) . Y: '5(z„w,.H(*),.) 

and times Tq"' < fj"' < . . . . We would like to define "^t — lini„_K>o however. 

In] 

the type processes Ku may not converge without some additional assumptions 
regarding the behavior of the migration processes {Z^, u g U}. 

Henceforth, we will also assume the following, which will ensure that for all 
n £ N and t > the expectation F[\{u £ U : Zu{s) £ En for some < s < ^}| ] is 
finite. 

Assumption 5.4. The sequence {En} of open sets can be chosen so that 
P"'We„ <t} <oo, for aU n e N and t > 0, 

where 

aA := inf{t > : Z{t) eA},Ae £. 



Remark 5.5. By our duality assumption, the measure m is stationary for Z. It 
follows easily that if for A e £ the condition P^{aA <i}<oo holds for some t > 0, 
then it holds for all t > 0. Furthermore, the condition P™{cryi < < oo for alH > 
is also equivalent to _P™[exp(— Act^i)] < oo for all (equivalent ly, s ome) A > 0. Using 
this equivalence, the question of whether or not Asssumption ^.4| is satisfied becomes 
a standard question in capacity theory. Under our duality assumption and the 
Assumption 5.1 that Z, Z are Hunt, Assumption 5.4 will certainly hold (with {£'„} 
any increasing sequence of relatively compact open sets such that 1J„ En — E) if the 
Lusin space E is locally compact and A-excessive functions for both semigroups {Pt\ 
and {Pt] are lower semi-continuous (see, for example, Remark 2.10 of |Get84]). In 
particular. Assumption p.4| holds if E is locally compact and Z and Z have strong 



Feller A-resolvent operators (see Exercise II. 2. 16 of [BG68|). Also, Assumption 5.4 
holds when Z and Z are Levy processes on R'' and m is Lebesgue measure (see 
Lemma 11.6 of [ |Bcr96t ). 



Fix i > and u E Un- Let 
(5.14) a„(u, t) := V sup{0 < s < t : Z„(s) 



Zi,{s) some v < u, v E Un}, 



and let /3„(u, t) be the corresponding value of u G Un, with Pniu, t) '■— u if a„(M, t) = 
0. Lemma 2.1 implies that /3„(u,t) is well-defined. In general, an(u,t) will not be 
one of the times {t],"'}, but we will have 



(5.15) 



= '^[r'(an(«,i)) = K2,^Aan{u,t)). 



Define /3^_„(s) := u for an{u,t) < s < t and /3^„(s) := Pn{u,t) for 
Cin{Pn{u,i),Cin{u,t)) < s < an{u,t). This definition extends iteratively to de- 
termine Pn „(s) on the interval < s <t with the property that 



(5.16) 



'it) 



is). 
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SO, in particular, Ku^\t) — f^^y Consequently, convergence of is equiva- 

lent to convergence of „ . 
For t > and u E U set 

(5.17) a{u,t) := V sup{0 < s < t : Z„(s) = Zy{s) some v < u, v Ell}. 

Let ?7 be a [0, l]-valued random variable that is cr (ZY)-measurable and takes values 
in the random set U (that is, C/ is a tT(i^)-measurable selection from U). By the 



duality assumption and the Hunt hypothesis Assumption 5.1 (cf. Proposition 15.7 
of |GS84[| ), t - a(U,t) has the same distribution as 



(5.18) inf{s > : Zu{s) = Zy{s) some v <U, v eU} At, 
where X^net/ \zu,u) is any Poisson random measure with mean measure 

(5.19) FxH^ I m{dz)P^{F)l{H) 

Je 

constructed from W using suitable further randomisation. Let f3{u,t) be the corre- 



sponding value ofv eU in ( p.l7D , with (3{u,t) := u if a{u,t) = 0. Define /^^(s) := u 
for a(it, t) < s < t and /?* (s) := /3(u, t) for a{(3{u, t),a{u, t)) < s < a{u, t). Extend- 
ing this definition iteratively, either we determine j3l^{s) on the interval < s <t 
and there are only finitely many levels in the range of /?* or there exists > 
such that limsjT^ Pui^) = 0- We show that this latter possibility cannot occur. 

Suppose that the latter possibility does occur. As above, let [/ be a criU)- 
measurable random variable taking values in lA. Define by analogy with T*, 
with the convention that := t \i there are only finitely many levels in the range 
oiPlj. Set 

(5.20) Z*,(r) :=limZ0*^(t_,,)(i-r'), Q<r<t-Tfj. 



Then by the strong Markov property, the duality assumption and Assumption 5.1 



{Zfj{r), < r < t — TIj} is a. cadlag Markov process with transition semigroup 
{Pr}- In particular, the range of this process is almost surely relatively compact 



and is contained in one of the En for n sufficiently large. Now, by Assumption 5.4, 
the cardinality of the set 

(5.21) {v e U : V < U, Zy{s) £ En some 0<s<t} 

is P-a.s. finite for all n E N. Consequently, P-a.s. there are indeed only finitely 
many levels in the range of and hence only finitely many levels in the range of 
/3* for all u gU. It follows that P-a.s. we have lim„^oo Pnu = ^^'^ hence 

(5.22) lim 4"l(t) = 4, (o)=:K„(i) 
for all u & U. 

If /3*^(s) = (3l^{s) for some < s < i, then /3^^(s') = Pl^is') for aU < s' < s. 
Moreover, if we define 

(5.23) Z*(r) = lim%(t_,,)(i-r'), < r < t, 

r' Ir 

for each u £ then conditional on U each Z* is a Markov process with transition 
semigroup {Pr}- In particular, {Z*, u g U} form a coalescing system of Markov 
processes, and for < r < t the equivalence relation defined by u ~ w if and 
only if /3*(t — r) = f3l{t — r) determines a partition {Ul,{r)} of the set of levels 
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U. For definiteness, assume that U is ordered U ~ {C/i, C/2, • ■ ■ } where the Ui are 
W-measurable random variables, and let Ul{r) be the equivalence class containing 
C/i, let Z^l(r) be the equivalence class containing the Ui with smallest index not 
contained in U{{r), etc. For each fc, .Z*(r) has the same value for all u € l^lir), 
which we denote by Zl{r). Then {{Z\,U{), . . . ) forms a coalescing Markov 

labelled partition oiU. 

Since the initial particle types {kJ^} are conditionally independent given {Z„, u G 
U} and U, and 



(5.24) 

for ui 
(5.25) 



^j,{Zui{0),dki) ■ ■ ■ ^i{Zu„iO),dkn)g{ki, . . . ,fc„), 

,Un eU, we have 

/ (Z„, (i ) , (t ) , . . . , (t) , K„„ (i ) ) I { Z„ (i ) } , W 
= HJ(Z„,(i),... ,Z„„(i)), 
where, in the notation of Section 0, 



(5.26) HJ(ei 



,e„) :-P 



/i(ZJ(i),dfcj)/(ei,fc^e(j), . . . ,e„,/c-ye(t)) 



for ei, . . . ,en £ E with ^ Cj, i ^ j. 
By and ( |5^ ) we have 



(5.27) 



X! f{Zm{t),Hm{t),-.- ,Zu„{t),Ku„{t) 

Ui, . . . ,Un 



f m®"(de)Ht/(e), 



for / a bounded measurabl e func tion on {E x if)". This identity gives a duality 
in the sense of (4.4.36) of EK86| between the discrete-particle, continuum-sites 
model and the corresponding coalescing Markov labelled partition process. 
Write 



(5.28) 

and 
(5.29) 



(Z„(t),A=„(t))- 



ueu 



Set — (7{Xl : s <t). The levels of the particles are independent of , so if 
/(e, k, u) satisfies /g^jg ^ m{de) l{du) supj.^^ |/(e, fc, u)\ < 00, then 



(5.30) P 



Bx ATx [0,1] 



ExKx[a,l] 



Xl{de X dk) l{du) /(e, k, u). 
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Moreover, if /(e, fc, u) satisfies /g^ [o i] ''^i'^^) K'^'^) i^^Pi^'^PkGK /(e, A;, u)) — 1) < oo, 
then 



(5.31) 



exp 



ExKx[0,l] 



exp I / Xl{de X dk) log / l{du) exp{f{e,k,u)) 

'exK "'[0,1] y 



5.3. Measure valued, continuum sites, stepping stone model. We empha- 



size that Assumptions 5.1 and 5.4 are still in force. Consider A > 1. We increase 



the "local density" of particles in the above construction by replacing m by Am 
and select the levels to be i.i.d. uniform on [0, A] rather than [0, 1]. Following 
the construction of ^E*^ and above, define 



(5.32) 

and 
(5.33) 



^' A 



(Z„(t),K„(t))- 



Note that k„ only depends on locations and types of particles at levels v < u 
and we can construct from 'I''^ simultaneously for all 1 < A' < A by taking 
to be the restriction of 'I'^ to the particles with levels in [0, A']; that is. 



(5.34) 



E 



t — "(Z„(t),K„(t),«)- 

ueu'',u<\' 



Consequently, we may carry out the obvious construction to build '5^ and 
with levels in [0, oo[. The initial locations and the levels are such that 



(5.35) 



(Z„(0),K„(0),u) 



is a Poisson random measure with mean measure given by 



(5.36) 



Ax B xC< 



m{dz) n{z,B)l{C), 



where / is now Lebesgue measure on [0, cx)[, and for 1 < A < oo each of the can 
now be defined via ( |5.32| ) with := {u e : u < A}. The analogue of ( ^.27D 
becomes 



(5.37) 



Ui, . . . ,Un £ 

Ut ^ Uj,i ^ j 



. ^«„(i),K„„(i)) 



A' 



/ m®"(de)HJ(e), 

JE" 



where Ht is defined as in ( 5.26| ). 
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Define = a{X^, ^-f - : s < t), and note that, as in ( |5.30| ) and ( |5.3l| ), we 



have 

(5.38) P 

and 
(5.39) 



ExKy.[0,\] 



Bx/f x[0,A] 



X^{de X dk)l{du) f{e,k,u) 



exp / d^^f 

W£;xA'x[0,A] , 



exp 



( / X^{de X dfc) Alog ( 1 + Y / 

\JexK \ J[0, 



A] 



l{du) (exp(/(e,fc,'u)) - 1) 



Suppose that /(e, k,u) — for u > Xq. Then for A > Ao, the random variables 
on the left of (5.38) and (5.39) do not depend on A. Since for fixed t > the a- 
fields are decreasing in A, the left sides of ( ^.3^ ) and ( 5.39] ) are positive, reverse 
martingales, and hence converge P-a.s. as A t oo. It follows that converges 
P-a.s. to a random measure X^ satisfying 



(5.40) 



and 



(5.41) 



BxKx [0,oo[ 
f 

ExKxlO,ool 



V 3 



X^{de X dk)l{du) f{e,k,u), 



exp 



exp 




d'frf 



V 3 



ExKx[Q,oo[ 



X^{de X dk)l{du) (exp /(e, fc, u) - 1) 



In particular, by (5.41), for each t >0, is a doubly stochastic Poisson process 
(that is, a Cox process) with random mean measure given by X^ (8) I. 
Dividing both sides of (5.37) by A" and letting A — > oo, we have 



(5.42) 



(BxK)" 



X^{dei X dki)---X^{den x dfc„) /(d, fci, . . . ,e„,A:„) 



/ m®"(de)HJ(e). 
Je" 



Note also that 
(5.43) 



X-(. X K) 



By Remark 3.1 we can regard the measure m{de)ii{e, dk) as an element of S (which 
we will also denote by /i) and the random measure X^ as a S-valued random 
variable. By Theorem 4.1, X^ has the same law as Xt under Q'' for each t>Q. 
In fact, it is not difficult to show that (X(°°, i > 0) is a Markov process with the 
same finite-dimensional distributions as X under Q'^. We stress, however, that we 
have only constructed X!^ as an almost sure limit for each fixed t > rather than 
as an almost sure limit in some space of cadlag paths. 
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6. Dissimilarity for the stepping-stone process X 

Suppose in this section that the reference measure m is finite. Without loss of 
generahty, we can take m to be a probabihty measure. 

Definition 6.1. Consider G S. For n = 2, 3, . . . define the rfi^-order dissimilar- 
ity of V to be the quantity 

/n 
TO®"(de) v{ei) ({k e A"" : % ^ h, for all j ^ (.}) . 
i=i 

Note that I > D2{v) > D^^iv) > ■ ■ ■ > 0. Write b{v) := sup{n : D,,{v) > 0} for 
the maximal dissimilarity of where we set sup0 = 1. 

As we remarked in the Introduction, it is possible, by exactly the same argument 
used in Proposition 5.1 of | Eva97|| , to show that if Z (and hence also Z) is a 



symmetric a-stable process on the circle T with index 1 < a < 2, then for fixed 
t > there Q'^-a.s. exists a random countable subset {ki, fc2, . . . } of the type- 
space K such that for Lebesgue almost all e G T the probability measure Xt{e) is a 
point mass at one of the ki. Indeed, under suitable hypotheses a similar argument 
should extend to certain other processes for which points are regular. It is clear 
that if Xf also has finite maximal dissimilarity Q'^-a.s., then the set {ki, fc2, . . . } 



is, in fact, finite Q'^-a.s. Theorem 6.4 below provides a sufficient condition for 
the maximal dissimilarity D{Xt) to be finite. We are able to verify this condition 
when Z (and hence also Z) is Brownian motion on T (see Corollary p.3| and the 



beginning of the proof of Theorem 10.2). We suspect that the condition is also true 



for Levy processes on T for which points are not essentially polar (see |EP98| for 
an indication that this might be so). 

Definition 6.2. Observe that if n' > n, then 

(6.2) ((z["^...,zKl),(7l"'l,...,7lr'l),?E) 

has the same distribution as (Z^"! , 7["1 , ^["1), where we write 

striction of the partition ^'"'(t) to N„. Consequently, on some probability space 
(rit°°l, jF[°°l,p[°°I) there is an -valued process Z, an N°°-valued process 7, and a 
process ^ taking values in the space of partitions of N such that, in an obvious nota- 
tion, {{Zi, ... , Z„), (71, . . . , 7n),^|N„) has the same distribution as (Zl"l, 7!"!, ^1"1). 

Remark 6.3. Recall the definition of Z^"!, Zl"! and ^["1 from Remark O Let 
Z["1I and ^t"!-'- be defined from Z^"! in a similar manner to Z^"! and , with 
the difference that when two coordinate processes of Zl"! collide, rather than the 
one with the higher index being killed, a colliding particle is killed at random 
independently of the past (with both possibilities equally likely). It is immediate 
from the strong Markov property that (Z["1,^["1I) has the same distribution as 
(Z["l,^["l) for all n e N. Consider t > and a bijection /3 : N ^ N and define 
^^^\t), a random partition of N, by i ~j(/3)(t) j if and only if f3~^{i) ^^(t) /3~^(j). 
Then ((Z^(i)(0))igN, ^^^^) has the same distribution as ((Zi(0))ieN, £,)■ particular, 
for each t > the random partition i^(i) is exchangeable in the sense of Kingman's 



definition of exchangeable random partitions (see Section 11 of |Ald85|). 
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Theorem 6.4. For any /i e S and t > Q, the maximal dissimilarity D(Xt) under 
is stochastically dominated by the number of blocks in the partition ^{t). In 
particular, if for some t > the partition ^{t) has finitely many blocks ¥^°°^-a.s., 
then D{Xt) < oo, Qf'-a.s. 

Proof. Fix a diffuse probability measure k on K (for example, k could be fair 
coin-tossing measure). Given another probability measure p on K, let de- 
note the push-forward of the measure p k on K x K hy the mapping (fc, h) i— > 
{ki, /ii, ^2, ^2, ^3, ■ • • ) from K X K into K. Let p^ denote the push-forward of p 
by the mapping fc i— > {ki, k^, k^, . . .) from K into K. Thus the operations p p^ 
and p I— *■ p^ are one-sided inverses of each other: we have (p~*)^ = p- Given 
/i £ S, define S S by p^{e) := fJ-{e)^ and p^{e) := p,{e)^ . Of course, the 

operations p^ and p i—>- p^ are also one-sided inverses of each other. Note 
for any /i S S that p~*{e) is diffuse for all e E E. 



Fix t>0. It is straighforward to check from the definition in Theorem 4.1 that 
the distribution of under Q'^ coincides with the distribution of Xt under Q^ . 
(This is, of course, what we expect from the stepping-stone model interpretation: 
a model that keeps track of the types for one trait should look the same as a model 
that keeps track of the types for two traits if we don't look at one of the traits.) 
Clearly, Dniv^) < Dn{i') for any n and S S, so D{v^) < D{i^). Consequently, 
D{Xt) under is stochastically dominated by D[Xt) under . 



We can use Remark 4.4 to compute multivariate moments of the form 



W [{An {Xt)Y' ■ ■ ■ {Dn, iXt)r] , G {2, 3, . . . }, a, e N, 1 < i < £, ^ e N, 

and discover that they are independent of ly within the class of G S with the 
property that i/(e) is diffuse for all e G E. Because < Dk{Xt) < 1 for all fc > 2, 
the multivariate moment problem for each of the vectors (Z3„j(Xt), . . . ,D„j(Xt)) 
is wcll-posed and hence the joint distribution of {D2{Xt), D^{Xt), . . . ) under Q'' is 
the same for all such v. Consequently, the distribution of D{Xt) under is also 
the same for all such v. In particular, if A G S is defined by A(e) := k for all e G i?, 
then the distributions of D{Xt) under and Q'^ are the same. 

Putting the above observations together, we see that it suffices to show that 
D{Xt) under is stochastically dominated by the number of blocks of £^{t). 

Let (Li)igN be an i.i.d. sequence of if -valued random variables (which we sup- 
pose are also defined on {Q^°°^ , , )) that is independent of Z with Li having 
distribution k. Put Li = ^^.^(t), so that Li — Lj if and only if i ~{(t) j, P[°°l-a.s. It 



follows from Remark 6.3 that the sequence ((Zi(0), Li))igN oi Ex if -valued random 
variables is exchangeable. 

Let Ai denote the point mass at {Zi{0), Li). By an extension of the standard 
reverse martingale proof of de Finetti's theorem, as n —^ oo the sequence of random 
probability measures ¥„ n^^ X]r=i converges p[°°l-a.s. in the weak topology 
to a random probability measure Y on E x K. Moreover, if we let T denote 
the permutation invariant cr-field corresponding to ((Zi(0), Li))igN (that is, X := 
fin <^{Y^n, Yn+i, •••}), then we have 

(6.3) p[°°l[0((Zi(O),Li),... ,(Z„(0),L„)) / 
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for any bounded Borel function (j) on {E x K)". (See the proof of Theorem 2.4 of 
|DK96 | for the details of this sort of argument.) 

We claim that Y has the same distribution as Xt under (recall from Remark 
3.1 that we can identify j/ € 5 with the probability measure 7n{de)v{e){dk) on 
E X K and that S-valued random variables become random probability measures 
on E X K when thought of in this way). If is a bounded Borel function on 
{E X K)" for some n e N, then, by (|0), 



dY 



(g)n 



.P[°°] [0((Zi(O),Li),... ,(Z„(0),L„))] 



jerH(t) 



K{dk,)cj, ((4"'(0), fc^H(,)), . . . , (4"1(0), fc^H(,))) 



A(z]"l(t))(dfc,)0((zi"l(O), 



, (4-1(0), A:^H(,)) 



Comparing this with the equivalent definition of {X, Q^) in Remark 4.3 shows that 
Y does indeed have the same distribution as X* under O'*'- 



(6.4) 



Finally, by (6.3) we have 

= 1 {34, . . . : if,., 1 < ^ < J < n} 

= 1 {i{t) has at least n blocks} . 



It is thus certainly the case that D{Xt) under is stochastically dominated by 



the number of blocks of ^(t). 



□ 



Remark 6.5. If Z and Z are both Hunt processes (that is, if Assumption holds), 
then the particle representation of Section ^ can be used to give a somewhat more 
direct proof of Theorem (note that Assumption 5^ holds because m is a prob- 
ability measure). We can sketch the proof as follows. The set of levels U°° in the 
construction of Section || is the set of points of a Poisson random measure on [0, oo[ 
with Lebesgue intensity and hence U°° is discrete. The dissimilarity i?„(X(°°) is 
just the conditional probability (conditioning on X!^) that the particles with the n 
lowest levels are all of different types. The argument that lead to (5.27) establishes 
that the total number of types exhibited by all particles is just the number of blocks 
in the corresponding coalescing Markov labelled partition. 

Remark 6.6. In the spirit of the previous remark, it is easy to see that D[Xt) is 
almost surely finite for all t > when Z (and hence also Z) is Brownian motion on 
the circle T and m is normalised Lebesgue measure. Once again, we just sketch the 
argument as a more quantitative result will follow from Corollary 9.3 below. Almost 
surely, there will exist two particles, say with the m"^ and v}^ lowest levels, m < n, 
such that by time t these two particles have collided and after the collision the first 
particle moved around the circle clockwise while the second particle moves around 
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anti-clockwise until they collided again. The total number of types exhibited by 
all particles at time t is then at most n — 1. 

7. Sample path continuity of the stepping-stone process X 

Our aim in this Section is to pre sent a sufficient condition for X to have con- 
tinuous sample paths (Theorem and use it to establish that if the migration 
Markov process is a Levy process or a "nice" diffusion, then X has continuous sam- 



ple paths (Corollary 7^ , Corollary |7.4| and Remark |7.5| ). The proof of Theorem ^2 
is postponed to the next section. We emphasize that we are no longer assuming 
that the reference measure m is finite. 

Definition 7.1. For e = (ei, 62) G E"^ with ei ^ 62, let T" := ini{t >Q: Z"^ {t) = 
Z^-^it)} denote the first time that Z''^ and Z''^ coUide. 

Theorem 7.2. Suppose there exists e > such that for all non-negative ■0 G 
L^{m) n L°°(to), 



limsup^-" / m®2(^e)^®2(g)p|ye < ^| ^ 
tio J 



UO 

Then X has continuous sample paths Q^-a.s. for all /i G S. 

Corollary 7.3. Suppose that Z is a Levy process on or the torus for some 
d G N, and m is Lebesgue measure. Then X has continuous sample paths Q^^-a.s. 
for all fi €z S. 



Proof. For d > 2 we have that T'^ = 00, P-a.s. for m®^-a.e. e, and so Theorem 7.2 



certainly gives the result. In fact, it follows from the remarks at the beginning of 



Section 5 in Eva97 | that X evolves deterministically and continuously in this case. 

Now consider the case where Z is R-valued. The T-valued case is similar and is 
left to the reader. 

Write {Z, P^) for the Levy process that is the symmetrisation of Z. That is, the 
distribution of Z starting at is the same as that of Z' — Z", where Z' , Z" are two 
independent copies of Z both started at 0. Put 

(7.1) T° mi{t > : Z{t) = 0}. 

Then for non-neg ative ?A e L'^{m) n L°°{m), 

m'^^{de)i;^\e)P{T'' < t} ^ f m®^{de)il;'^'^{e)P^^-^^{f° < t} 



(7.2) 



m{dx)^l;{x)P''{T' < t}, 



where ip{x) := J m{dy)ip{x + y)'ip{y) G (m) L°° (m) and we are, of course, using 
the shift invariance of m. 

For a > write C", U°' and e" for the a-capacity, a-resolvent and a-energy 



corresponding to Z (see Sections L2, ILS and n.4 of |Ber96| for definitions). Because 



Z is symmetric, these coincide with the corresponding dual objects. Write ^' for the 



characteristic exponent of Z (see Section LI of [Ber96|). Note that ^E* is real- valued 
and non-negative. 
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Using the convention ^ = 0, we have from Theorems II. 7 and 11.13 of | Bcr96| 
that 

(7.3) J m(d.)^-(x)P^[exp(-aTO)] = C"({0})C7"V5(0) < < 

By Proposition 1.2 of per96t , 

(7.4) ^{z) < cz^, \z\ > 1, 
for a suitable constant c, and so, for a > 1, 

(7.5) ae"({0}) = — / — dz > — dz > c'a^ 

for a suitable constant c'. 

Use the inequality l[o^t] i^) < e~°'^ + 1 — e^"* < e^""^ + at and take a = 3 to 
get, for t < 1, 

(7.6) / m{dx)i>{x)P^f'> <t}< c"||V;||ooa"^ + Miat < c*(||^||i + Moo)tK 



for suitable constants c" and c*. Now apply Theorem 7.2 with e = 1/3. □ 

Corollary 7.4. Let d be a metric inducing the topology of the Lusin space E . Write 
B{x, r) := {y (z E : d{x, y) < r} for the closed ball of radius r > centred at x d E 
and := mi{t > : Z{t) ^ B{Z{0),r)} for the time taken by Z to first travel 
distance r from its starting point. Suppose that there are constants a, /3, 7 > such 
that 



and 



limsupr " sup m(_B(a;, r)) < 00 



limsupr supP^{S"' < r} < 00. 

rj,0 x£E 



Then X has continuous sample paths Q^-a.s. for all fj, ^. 
Proof. For non-negative ijj G L^{m) n L°°{m) and 5 > we have 

m®2(^e)V'®^(e)P{T'' < t} 
- J m®2(de)?A®2(e)i{d(ei,e2) < S}F{T'' < t} 
+ [ m®'^{de)iP®'^{e)l{d{ei,e2) > S}P{T'' < t} 



(7.7) 



< ll^llillV'lloo sup m{B{x, 6)) + 2\\ij\\l sup ^"{5*/^ < 

x£E xeE 



Take S = t^^ to get that the hypothesis of Theorem ^2 holds with e = (a/3) A 7. □ 

Remark 7.5. The above result can be applied to the case where Z is a regular 
diffusion on M in natural scale. In this case m is the speed measure and Z ~ Z . If 
m{dx) = a(x)dx with a bounded away from and 00 and d is the usual Euclidean 
metric on M, then it is not difhcult to see that the conditions of the corollary hold 
for a = 1, /3 < 1/2, and any 7 > 0. We leave the details to the reader. 
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8. Proof of Theorem 7. 



The proof of Theorem 7.2 will involve checking Kolmogorov's criterion for the 



sample path continuity of real-valued processes of the form {I{Xt; 4>))t>o for suit- 



able (j) G C{K), where /(•;•) is defined in Definition 3.2. The proof will be via 
several lemmas. 

Remark 8.1. In performing the necessary moment computations for Kolmogorov's 
criterion we will need to consider the various orders in which particles can coalesce 
in the coalescing system that define these moments and estimate the contribution 
of each possible sequence of collisions. We will repeatedly use the fact that if, for 
fixed i ^ j, we "swap" Z^^{t) and Z'^^{t) immediately after a stopping time S for 
Z*^ at which Z'^'(S') = Z^^iJS) to form a new process Z'', then Tf^ has the same 
distribution as Z*^. More precisely, if we define 



(8.1) Ztit) := 



(8.2) ^{t) 



(t), for t < S, 

ZJ{t), for t > S, 

'Zj{t), for t<S, 

Zf{t), foi t>S, 



and 

(8.3) Zf,:=Zl h^{t,j}, 

then, by the strong Markov property, Z'^ has the same distribution as Z®. 

Definition 8.2. For n',n" e N consider e' e E"' and e" e E"^" with 
e[, . . . , e^,, e", . . . e^„ distinct. Define a process taking values in the col- 

lection of finite sequences of two element subsets of {e^, ... , ej^,, e", . . . e'^„} and 

stopping times = T^"^ < ''^ < ... as follows. For e- € {e[, . . . , e'^,}, write 
S'^y M{t > : Z^'it) = Z"'^ (t) for some j ^ i 

or Z"'- {t) = Z^'i^ (t) for some k such that zf (s) 7^ f for all s <t}. 
For e'/e{ei',... ,e::„} write 



inf{t > : Z"'^ (t) = Z^'^ (t) for some j 
or Z^'^it) ^ Z'''^{t) for some k ^ i such that zf {s) ^ ffor aUs < t}. 

Loosely put, we are thinking of the particles starting at coordinates of e' as evolving 
freely without coalescence whereas the particles starting at coordinates of e" are 
undergoing coalescence among themselves. Moreover, St ''^ (if it is finite) is the 
first time that the particle starting at e (where e is either a coordinate of e' or e") 
collides with a "living" particle starting at one of the other coordinates. 

Let < n' + n" denote the cardinality of the random set of time points 

(8.4) : e G {e'l, . . . , e^,, e'/, . . . e'„'„} and S^'l''" < 00} 

and, if R^'^^" > 0, write T^^^ < . . . < T^J'*|e" fo'' an ordered listing of this set. 
Put Tn''*'" and Tf''*'" := 00 for £ > R'''^''" . 
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Set $^'l*'"(ro''''''") := 0. For 1 < fc < write 
(8.5) {xfe,yfc}C {e'„... ,e;„e'/,... ,e;'„}, 

for the P-a.s. unique unordered pair such that Z^''{T^ ) — Zy''{T^ ). By 
definition, at most one of Xk and yk belong to {xi,yi, . . . , Xfc-i, Put 
^e'|e"^ye'|e"^ := ({xi,yi},... ,{xfc,yfc}). Complete the definition of by 
setting $''>"(i) := $''>"(r/'''"), where ^ > is such that T/''''" < t < T/|'f' . 

Definition 8.3. Given e e n E N, with ei, . . . , e„ distinct, define a process 
^'^ taking values in the collection of finite sequences of two element subsets of 
{ei, . . . , e„} and stopping times — Tq < Tf < • ■ • by (with a slight abuse) re- 
using the definitions of and ''^ with e' = e and e" the null vector. That 
is, all particles evolve freely with none of them killed due to coalescence. Note that 



if n = 2, then Tf = T*^, where T"^ is the first collision time from Definition 7.1 



Notation 8.4. Given e' e E"- and e" e i?" , write e' : e" for the concatenation of 
these two vectors. That is, e' : e' := {e[, ... , e'„,,e'{, . . . e'^„) £ £;"'+"". 

Notation 8.5. For x G M write \ x\ for the greatest integer less than or equal to x. 

Lemma 8.6. For non-negative ip G L^{m) H L°°(m), t > 0, and q,n',n" Q N we 
have 



m®"' (E) m®""(rfe' (g) rfe")V'®"'(e')'/'®""(e")P{r,'''l''" < t} 
(8.6) < / TO®"' (g) m®"" ide' (g de") V^"' (eOV-®"" (e")P{T/^^" < t} 



< c{n' + n", q, ^P) (^j m®2(rfe)V®^(e)P{r^ < t}^ 
for some constant c{n' + n" , q, ip) that depends only on n' + n" , q and %p. 

Proof. By definition, S^. < Se for any e G {e[, . . . , ej^,, e", . . . e"„}, and so 
jie :e |e ^j^j therefore suffices to show that for g, n G N 

J m'^"{de)iP®'\e)P{T^ < t} 

(8.7) ^ L§J 

< c(n,9,V) ( / m®2(^e)^®2(g^p|j.e < ^1 



We begin with some notation. For any sequence of pairs 

H = {{xi,yi}, ... , {xi,yi}), with Xi.y^ G E, and x^ 7^ for I < i < \H\ ^, 
and t > 0, define an event 

(8.8) :— {T^"^^'^^^ < ^(^2,^2) < _ . < j^ixe.Vi) < 

It is easy to see that Af D Af^ for any subsequence H' of Put H := 

ULd^i^y^} ^ ^- For z G 17, define ^(z, iJ) = {1 < i < : z G {x^,y^}} to 
be the set of indices of the pairs in which z appears. 
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^ Now fix G = {{xi,yi},... ,{xg,yq}), with Xi,yi e E, ^ y^, 1 < i < q, and 

G ULii^i' Vi} ^ {"^i' ■ ■ ■ ' S"}- We wish to estimate f{Af). 

Let . . . , i/i) be the subsequence of (1,2,... , q) obtained by hsting the ele- 
ments of {max(.(2;, G), 2; £ G} in increasing order. Define a subsequence G* of G 

by 

(8.9) G* := ({a;,i,y^J, . . . , {x;,^ , j/ij) =: ({xo,!, ?/o,i}, ■ • • , {a;o,|G.hyo,|G.|}) 

(the reason for the alternative indexing will become clear as we proceed) . Note that 
|G| = |G,| because for ah z e G, z e {a;,nax 1(2,0), ?/max 1(2,0} ^ '^*- definition, 
for 1 < j < |G*| the inequalities maxi(a;ij , G*) > j and max(.(?/i^-,G*) > j hold, 
and at least one of these inequalities is an equality. In other words, 

(8.10) min{maxt(xoj,G*),maxi(j/o,j,G*)} = j, l<j<|G*|. 

Without loss of generality we can assume that zi = max ^.(2;^^ , G) < maxi(yij,G). 
Then Xi^ ^ {xr^yr} for ii < r < q, and, a fortiori, Xi^ ^ {xi^^yij,} for 1 < p < h. 
Hence, i(xo,i, G*) = t(a;ii j G**) = {1} and we are now in one of the following three 
cases: 

Case I: |(-(yo,i, G*)| = 1. Let Gi be the subsequence of G* obtained by deleting 
{xo,i,yo,i}. Then Gi n {a;o,i, yo,i} = 0, 

P(Af ) < P(ylf *) < P [{T^^o.uvo.i) < ^j. n ^Gi^ 
and |Gi| ^ |G;| - 2= |G| - 2. 

Case II: |i(?/o,i, G*)| — 2. Write t(yo,i)G*) = {l,j2}- Define a subsequence Gi of 
G* by deleting {xo,i,?;o,i} and {2:0,^2,2/0,^2}) from G*. Then Gi n {xo,!, 2/0,1} = 0, 

P(Af ) < P(Af *) < P [{T^^o.uyo.i) < n ^Gi^ 
= P{T(^»-1'2«'-i) < i}P(4^i), 

and |Gi| > |G,|-3= |G| - 3. 

Case III: |t(yo,i,G*)| > 2. Write i(yo,i,G*) = {1,^2,- ■• ,jp} where 1 < j2 < 
... < jp and j/0,1 = Va,]2 = ■■■ = Vo,]^- Then maxi(a;o,j2 , G*) = j2 because 
maxi(yo,i2, G*) = jp > j2- Let 

(8-13) G*, := i{xQ.i,yo^i}, . . . ,{a;o,|G.h2/o,|G.|})> 

where 



•14) y'o. 



yo,j, if J < J2 or 2/o,i 7^ yo,i, 
Xo,]2 , if J > J2 and ?/o,j = 2/o,i- 



We then have G* = G** and by, switching Z^'^^a and at time T^^°-^2'V0:32) 

in the manner described in Remark 8.1, we also have P(Af*) = ¥{Af"). More- 



over, t(xo,i,G,,) = 1 because a;o,i ^ U?=,2{^o,i, 2/o,j} = UL,2't^oj'2^o )}• 
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|t(yo,ij G**)! = 2 and we are in Case II with G, replaced by G»». From the dis- 
cussion in Case II we know that there exists a subsequence Gi of G»* such that 
Gin{a;oa,2/oa} = 0, 



(8.15) 

and |Gi| > |G, 



|G, 



|G|-3. 



The reduction procedure that transformed G into Gi can be repeated at least 
([|G|/3J — l)"*" more times. That is, for < £ < [|G|/3J, there exist sequences 
Ge = {{xe^i,yes}, ■ ■ ■ ,{xe,\Ge\^ye^Gi\} such that Go = G and for < ^ < [|G|/3J-1 

(8.16) P{Af'-) < P (^{T^^^'i'^'^'i) <t}n '+^) ^ p{r(^i,i.m,i) < t}P{Af'+'), 

with Gi+i n {xe,i,yis} = 0, G,+i C G^, and |G,+i| > \Gi\ - 3. 
It follows that 

(8.17) P(Af ) < p{r(^o,i,yo.i) < t}¥{Af^) <...< Yl P{r(^'-i^2'''i) < t}, 

where the sets {a;i,i, i = 0, 1, . . . , [|G|/3J — 1, are pairwise disjoint. 

Write Gf for the set of possible values of ^^{T^) on the event < t. Note that 



|G| > q for any G G C and so the rightmost product in (8.17) has at least lq/i\ 
terms. Therefore, if we let c(n, q, -0) denote a constant that only depends on n, q, 
(but not t), we have from (8.17) that 



m'^"{de)tl;'^"{e)P{T^ < t} 

= J m^"(de)V®"(e) ^ P{$*'(T,'=) = G, < t} 



(8.18) 



GeGf 



< / m®"(de)i/.®"(e) ^ 
GeGf 



< c(n,g, V) / m^''(de)V®^(e)P{r'= < i} 



LfJ 



□ 



(8.19) 



Notation 8.7. Given /i G ^, define /it G S, t > 0, by 

/(/it;,/.) := / m{de)P J fi{Z''{t)){dk)cj>{e){k) 

m{de)F J /i(e)(dfc)0(Z"(i))(fc) 
so that/(/it;0)=Q''[/(Xt;0)]. 
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Lemma 8.8. Consider /i e S, i > 0, g G N, and <j) — ^ ® x with non-negative 
■0 e L^{m) n L°°{m) and non-negative x ^ ^i^)- There exists a constant c{(l),q) 
that only depends on <f> and q ( and not fi or t) such that 



>^ {liXt; cf,) - /(Mt; < c{4>, q) [ I m«-^(de)V'«^(e)P{T- < t} 



LfJ 



Proof. Given n',n" G N and vectors f G and f" G E" with f[,...J'^,, 
/{',.•■ distinct, write 



.20) 



L{zf''{mdK,)®6^ii{zf^"{t)){dk', 



i' = l 



i" = l 



V®"' (h' h' \ V®"" ( h" h" \ 



Then, by definition, 



2q 



2q 



Therefore, by Lemma 8.6, it suffices to show that 



.21) 



m®2«-*(de') ® m'^\de")il)"^'"i-\e')ilj"^\e") 



i=0 
X 

= 



[L^'\^"{t)i[T;'\^" <t<T;li] 



for < p < g — 1. 

For e' G E'^'?-* and e" G with e'^, . . . , Cj^^j, e", . . . , e" distinct, write 

Sji'' 1 < j < 2q — i, < ft. < i, for the coUection of subsets of 
{e'l, . . . , e2^_j, e'/, . . . , e"} with exactly j elements from {e'l, . . . , e'a^^j} and exactly 
/i elements from {e", . . . , e"}. Put 



.22) 



2(7 - i 



, 0<j<2q-i,0<h<i. 
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It is clear by construction that, recalling the transformation H is- H from the proof 
of Lemma B.6, 



= Z1 / m®^'^-\de')(g,m^%de")iP^^'i-'{e')iP®%e") 



X 



Note that |<i)'''l®" (Tp )| < 2p, and so a necessary condition on j, h for a summand 
in the last term to be non-zero is that j + h < 2p < 2q. 

For fixed e' and e" write e* := (e^,... ,e2q_i,e'/^^i, . . . ,e'') and e** := 
(e", . . . , e'^). Observe that 

{t;>" < t < T;'lf, $^(T;>") = {e'l, . . . , e;, e'/, . . . , e;,'}} 

(8.23) ^ „ J 

= {t;I- <i<r;Jf , <i>e*l«"(T; I'' ) = {el,... ,e;,e'/,... ,e^}}. 

Moreover, on this event the partition (t) is the restriction of the partition ^® (t) 
to N/i, and hence L"^ I'' (t) = U' (t) on this event. Therefore, the quantity 

(8.24) 

m®2'?-'(de') ® m^\de")i;^^''-'\e')^^\e") 



i">"(t)i {r;>" < < < T;}f, (t/I'=") = {e'l. 



/ ff ff 



>} 



does not vary as i ranges from h to 2q — j. 

The proof is complete once we note that for fixed h, j with h < 2q — j we have 



.25) 



(2g)!(-l)'' (2g-J-/i)! 



(2g-j-/i)!j!/i! 



i—h 



{z^hy.i2q-i-jy. 



(2g)!(-l) 



(29-j-/i)!j!W 



(1 _ l)29-j-/« 



0. 



□ 



Completion of the Proof of Theorem 7.2. Because X (as a Hunt process) has 
cadlag paths Q^-a.s. for all e S, it suffices to show that I{X.\(f)) has continuous 
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sample paths Q^-a.s. for all G S and all 4> belonging to some countable subset 
of L^{m,C{K)) that is separating for S. Moreover, because I{X.;(j)) already has 
cadlag paths, verifying Kolmogorov's criterion establishes that these paths are, in 
fact, Q'^-a.s. continuous. That is, verifying Kolmogorov's criterion does more than 
just establish the existence of a continuous version of X, it establishes that the 
version we already have is continuous. 

Let {C/"}a>o denote the resolvent corresponding to the semigroup {Pt}t>o- Sup- 
pose that iS is a countable collection of bounded, m-integrable, continuous, non- 
negative functions on E with dense linear spa n in L^ {m) (such a collection can 
be seen to exist by combining Lemma A.l of |Eva97| with Proposition 3.4.2 of 
EK86t ). Note that if 6* e 5, then aU°'0 converges to pointwise as a — > oo. 
Also, / m{dx)aU'^9{x) — J m{dx)9{x) < oo by the duality hypothesis for the 
pair Z, Z . By a standard extension of Lebesgue's do minate d convergence theorem 
(see, for example. Proposition 18 in Chapter 11 of Roy68 |), if g G L°°{m), then 
lim^^oo / m(dx)aU°'9{x)g(x) = J m{dx)9{x)g{x). 

Write V := {U"9 -.0 € S, a rational} C L^{m) n L°°{m). It follows easily from 
what we have just observed that if C is a countable dense subset of {x & C{K) : 
X > 0}, then the countable collection of functions of the form ij} ® Xi with ip E T) 
and X G (5, is separating for S. 

Fix ip € T) (with = U°'0 for 6* e 5 and a rational), x G C, and g G N such 
that [(//Sje > 1, where e > is as in the statement of the theorem. In order to 
show that I{X.\ ip'S)x) Q^-a.s. continuous sample paths for all /i, it sufhces by 
the Markov property of X and Kolmogorov's continuity criterion to show for some 
constants c and S which depend only on ?/;, X) 9 that 



.26) 



{i{Xt;^(g)x)- Hp;^(^x)} 



2q 



< ct 



l+S 



for alH > and /i G S. This, however, follows from Lemma S.8 and the observation 
that 



|/(/it;V'«'x) - ^(m;V'«)x)I < I m{dx) Pt^p{x) - ipix) 



m{dx) 



< 2a"^(l-e""*) / m{dx)0{x) < 2t m{dx)0{x) 



.27) 



where we have used the consequence of the duality hypothesis on Z, 
^ m{dx)aPt0{x) = ^ ra(dx)0{x). 



Z that 



9. Coalescing and annihilating circular Brownian motions 

In this section we develop a duality relationship between systems of coalescing 
Brownian motions on T, the circle of circumference 27r, and systems of annihilating 
Brownian motions on T (Proposition pj] ). This relation will be used in Section [lO| 
to investigate the properties of the stepping-stone model X when the migration 
process is Brownian motion on T. It will also be used in Section |l^ to study the 
random tree associated with infinitely many coalescing Brownian motions on T. We 
mention in passing that coalescing Brownian motion has recently become a topic 
of renewed interest (see, for example, [ TW97 | and [PsiOS ). 
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For the rest of this paper, Z (and hence Z') will be standard Brownian motion 
on T, and m will be normalised Lebesguc measure on T. 

Given a finite non-empty set ACT, enumerate A as {ei,... ,e„}, put e := 
(ei, . . . , e„), and define a process W-^, the set-valued coalescing circular Brownian 
motion, taking values in the collection of non-empty finite subsets of T by 

(9.1) W^{t) :={Z^e(,)(0,... ,^^e(,)(i)}={^^e(,)(i),... , Z^^t^t)} , t > 0. 

Equivalently, (t) is the set of labels of the coalescing Markov labelled partition 
process ('^{t). Of course, difi^erent enumerations of A lead to different processes, 
but all these processes will have the same distribution. In words, describes 
the evolution of a finite set of indistinguishable Brownian particles with the feature 
that particles evolve independently between collisions but when two particles collide 
they coalesce into a single particle. 

Write O for the collection of open subsets of T that are either empty or consist 
of a finite union of open intervals with distinct end-points. Given B & O, define on 
some probability space (S, Q, Q) an O-valued process , the annihilating circular 
Brownian motion as follows. The end-points of the constituent intervals execute 
independent Brownian motions on T until they collide, at which point they anni- 
hilate each other. If the two colliding end-points are from different intervals, then 
those two intervals merge into one interval. If the two colliding end-points are from 
the same interval, then that interval vanishes (unless the interval was arbitrarily 
close to T just before the collision, in which case the process takes the value T). 
The process is stopped when it hits the empty set or T. 

We have the following duality relation between and . 

Proposition 9.1. For all finite, non-empty subsets ACT, all sets B £ O, and 
all t>0, 

f{W^{t) C 5} = Q{A C V^{t)}. 

Proof. For G N, let := {0,1, . . . N - 1} denote the integers modulo N. Let 

:= |, . . . , } denote the half-integers modulo N. A non-empty subset 
D oi Zn can be (uniquely) decomposed into "intervals": an interval of D is an 
equivalence class for the equivalence relation on the points of D defined by x ~ y if 
and only if x = y, {x, a; + 1, . . . ,y—l,y} C D, or {y,y + l,... , a; — 1, a;} C D (with 
all arithmetic modulo N). Any interval other than Zjv itself has an associated pair 

of (distinct) "end-points" in Z^: if the interval is {a, a + 1, . . . ,6—1,6}, then the 
corresponding end-points are a — ^ and b+-^ (with all arithmetic modulo N). Note 
that the end-points of different intervals of D are distinct. 

For C C Zjv, let Wf^ be a process on some probability space (n',JF',P') taking 
values in the collection of non-empty subsets of Zjv that is defined in the same 
manner as W"^, with Brownian motion on T replaced by simple, symmetric (con- 
tinuous time) random walk on Zat (that is, by the continuous time Markov chain 
on Zat that only makes jumps from a; to a; -|- 1 or a; to a; — 1 at a common rate A > 
for all X e Zjv). For D C Zn, let be a process taking values in the collection 
of subsets of Zat that is defined on some probability space (S', Q' , Q') in the same 
manner as , with Brownian motion on T replaced by simple, symmetric (contin- 

uous time) random walk on Z^ (with the same jump rate A as in the definition of 

W^). That is, end-points of intervals evolve as annihilating random walks on Z^. 
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The proposition will follow by a straightforward weak limit argument if we can 
show the following duality relationship between the coalescing "circular" random 
walk and the annihilating "circular" random walk : 

(9.2) r{W^{t) CD}^ Q'{C C V^{t)} 

for all non-empty subsets of C C Z^v, all subsets of D C Z^r, and alH > 0. 

It is simple, but somewhat tedious, to establish ( |9.2[ ) by a generator calculation 
using the usual generator criterion for duality (see, for example, Corollary 4.4.13 
of pK86|| ). However, as Tom Liggett pointed out to us, there is an easier route. A 



little thought shows that is nothing other than the (simple, symmetric) voter 
model on Z^v- The analogous relationship between the annihilating random walk 



and the voter model on Z due to |Sch76] is usually called the border equation (see 



Section 2 of |BG8C| for a discussion and further references). The relationship ( p^ ) is 
then just the analogue of the usual duality between the voter model and coalescing 
random walk on Z and it can be established in a similar manner by Harris's graphical 



method (again see Section 2 of |BG80| for a discussion and references). □ 



Remark 9.2. We have been unable to find an explicit reference to Proposition 3.1 
or its analogue for Brownian motion on R. However, if, in the R-valued analogue, 
one considers a limit where A approaches a dense subset of R, so that the system 
of coalescing Brownian motions converges to a coalescing Brownian flow, then the 
analogous result for the limiting flow can be found on plS of [ Arr79| . 



Recall Z and 7 from Definition 3.2. Define set -valued processes W^"\ n & N, 
and W by 

(9.3) W["l(t) :={Z^,(4)(i),... ,Z^„(i)(0}CT, t>0, 
and 

(9.4) W{t) := {Z^,^t){t), Z^,^t){t), . . . } C T, i > 0. 

Thus, W^^Ht) C W^^^t) C ... and U„eN^'"'(i) = W{t). Recah that (l¥["l(i)t>o 
has the same law as ({z'"' (t),... , z'?' {t)})t>o- Put N(t) := \W(t)\, the 

7I (t) Tr^ '(*) 

cardinality of the random set W{t). Note that N{t) is also the number of blocks 
in the partition ^{t), which is in turn the cardinality of the random set T{t). It is 
clear that p[°°l-a.s. N{t) is a non-increasing, right -continuous function of t and if 
N{to) < 00 for some to > 0, then N{t) - N{t~) is either or -1 for all t > to- By 
the following corollary, N{t) < 00, p[°°l-a.s., for aU t > 0. 



Corollary 9.3. For t>0, 



and 



[A^(t)] = l + 2^exp(-(^)"'t) <(X3 



limtlpM \N(t)] = 2V^. 

no 



Proof. Note that if _B is a single open interval (and hence for alH > the set (t) 
is either an interval or empty) and we let L{t) denote the length of V^{t), then L 
is a Brownian motion on [0, 27r] with VarL(t) = 2t that is stopped at the first time 
it hits {0,27r}. 
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(9.5) 



= 1 
= 1 



Now, for M g N and < i < M ~ I we have from the translation invariance of 
Z and Proposition 9.1 that 

pM n [2tti/M,2tt{i + 1)/Af] 7^ 0} 

C]0,27r(M-- 1)/Af[} 

where we take the annihilating process 

y]o,2^(M-i)/M[ defined on the same 

probability space (f]M, jcM^pM-j 

process Z that was used to construct 
and W, and we further take the processes y]o,27r(M-i)/j\/[ ^^^^ 2 to be inde- 
pendent. Thus, 

pM n [2Tri/M, 2TT{i + 1)/M] 7^ 0} 

(9.6) = 1 - P[°°l |ylO,2^(M-l)/Af[(^) ^ 

= 1 - p|f < 2<, B{f) = 27r I 5(0) = 2tt{M ~ 1)/m}, 

where _B is a standard one-dimensional Brownian motion on some probability space 
(f),j£",P) and f = inf{s ^ : B {s) e {0,27r}}. 
By Theorem 4.1.1 of |Kni81| we have 



[\W{t)\] 
hm p[°°l 

M^oo 



'M-1 



1 {M^(<) n [2ni/M, 2n{i + 1)/M] ^ 0} 



lim M 

M^oo 



(^1 - P |f < 2t, B{f ) = 27r I 5(0) = 27r(M - 1)/m|) 



2 °° (—1)" 
1 — lim M— sin | nn 



I'M - 1 
\ M 



exp 



= 1 + 2 ^ exp 

n=l 



where 
(9.7) 



0(u) 



exp(— vrn^w) 



n— — 00 



is the Jacobi theta function (we refer the reader to | BPY98[ for a survey of many of 
the other probabilistic interpretations of the theta function). The proof is completed 
by recalling that 9 satisfies the functional equation d{u) = u~2 9[u~^) and noting 
that lim„_^oo ^(u) = 1- D 

We conjecture that t^N{t) 2^/tt as t | 0, p[°°l-a.s. However, we are only able 
to prove the following weaker result, which will be used in Section |l^. The proof 
will be given at the end of this section after some preliminaries. 
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Proposition 9.4. With p[°°l -probability one, 

< liminf t'iV(t) < limsupi3 7V(t) < oo. 

no t|o 

For t > the random partition ^{t) is, by Remark 3.3 and Corollary 6^3, ex- 
mgeable with a finite number of blocks. Let 1 = x\ < x\ < ... < x^,^-. be 



chan; 



the list in increasing order of the minimal elements of the blocks of ^(t) (that is, 
a list in increasing order of the elements of the set r(f)). Results of Kingman (see 
Section 11 of [Ald85] for a unified account) and the fact that ^ evolves by pairwise 
coalescence of blocks give that p[°°l-a.s. for alH > the asymptotic frequencies 

(9.8) F,{t) = lim n-^\{3 e N„ : j ^^(t) 

exist for 1 < i < N{t) and Fi(i) H h FN(t){t) = 1. 

Lemma 9.5. With ¥^°°^ -probability one, 

N(t) 



limi-5 FM"^ 

no ^ ^ ' 



r3/2 ■ 



Proof. Put M{t > : Z,{t) = Zj{t)} for i ^ j. Observe that 



N{t) 



^ n n 

1™ ~EEi 0' ^} 



P[°°1{1 2} 



From Theorem 4.1.1 of [KniSlj we have 
pM{Ti2 < t} 

= — 1 - - > si" ' ^ ' 



1 



= -y 



l(2n-l)2 



1 — cxp — 



2n- 1 



2n~ 1 

2 



exp 



2n- 1 



271 - 1 



ds 



2 /"* 1 f °° s °° 

= 2 5: ^-pHi)- E -ph 

^ n— — oo n— — oo 



ds 



where 6 is again the Jacobi theta function defined in (9.7). By the properties of 9 
recalled after (9.7), 



(9.9) 



lim t 2 1 



Nit) 



= limi-2pM{7^ = 
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Now 



p[oo] 



/Nit) 



^ n n n n 

ii = l 12 = 1 i3 = l i4=l 

~^(t) 2, 3 ~5(t) 4}, 



and so 



(9.10) 



Var 



{1 2, 3 4} - pl°°]{Ti2 < T23 < 0- 



Observe that 



'"^{Tia < i, T34 < t, Ti3 > i, Ti4 > t, T23 > t, T24 > *} 

< PM{1 ~«(t) 2, 3 4, {{1,2,3,4}} ^ ^W(i)} 

< p[°°l{Ti2 < t, T34 < t} 



and 



P[°°l{ri2 < T34 <t}~ P[°°1{T12 < t, r34 < t, Ti3 > i, Ti4 > T23 > t, T24 > t} 

< ^ ^ P'~i{ri2 < t, T34 < t, Tij < t}. 

i=l,2 j=3,4 

Thus 



(9.11) 



Var Pii^f ^ 2 ~«(t) 3 4} 



i=l 



+ ^'"^^12 < i, 734 < t, Tij < t}. 

i=l,2 j=3,4 



Put := \Z^{0) - Zj{0)\. We have 



4} 



{T12 < t, Ti3 A T23 < t, Ti4 A r24 A T34 < t} 
{T12 < t, Ti3 A T23 < t, Tu A T24 A T34 < t] 



(9.12) 



\ {D12 < t-^, {Di3 A i^23) < t^, {Di4 A ^24 A D34) < t^} 

+ pl~l{Di2 < ti, {Di3 A D23) < ti , {Dii A D2i A D34) <ti} 
< Yl ^ ^ij >th + P'°°> I max Dij <3tH, 

l<i<j<4 11<«<J<4 J 

where we have appealed to the triangle inequality in the last step. Because | < 5, 
an application of the reflection principle and Brownian scaling certainly gives that 
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the probability p[°°l{Ty < t, D,j > ii} is o(P) as i j for any a > 0. Moreover, by 
the translation invariance of m (the common distribution of the Zi(Q)), the second 



term in the rightmost member of ( |9.12 ) is at most 

P[°°^{|^2(0) - Zi(0)| < 3tt, 1^3(0) - Zi(0)| < Zti, |Z4(0) - Zi(0)| < iti} 
= P[~1{|Z2(0)| < 3<t, 1^3(0)1 < iti, |Z4(0)| < iti} 

6 

for a suitable constant c when t is sufficiently small. Therefore, 
PM{1^,(,)2^,(,)3~,(,) 4} 

(9.13) = p[°°l{{Ti2 < i, ri3 A T23 < Ti4 A A T34 < t) 
= 0(tt), as no. 

A similar argument establishes that 

(9.14) p[°°l{ri2 < i, r34 < < i} 0(tS), as no, 

for i = 1, 2 and j = 3, 4. 

Substituting (|]l|) and ( |9.14D into ( p.ll| ) gives 



(9.15) Var ^F,;(t)2 =0(if), as n 0. 

This establishes the desired result when combined with the expectation calcula- 
tion ( |9.9| ), Chebyshev's inequality, a standard Borel-Cantelli argument, and the 
monotonicity of Fi(t)^. □ 

We may suppose that on our probability space (f^M, jrM^pM) there IS a se- 
quence Bi, B2, ... of i.i.d. one-dimensional standard Brownian motions with initial 
distribution the uniform distribution on [0, 27r] and that Zi is defined by setting 
Zi{t) to be the image of Bi{t) under the usual homomorphism from M onto T. For 
n e N and < j < 2" — 1, let < /^'"' < ... be a list in increasing order 
of the set of indices {i e N : Bi{0) e [27rj/2", 27r(j + l)/2"[}. Put B"'^ := Bj,^., 

and Z"'-* := Zjn.j. Thus (_B"'-')igN is an i.i.d. sequence of standard M-valued 

Brownian motions and (Z"'-')igN is an i.i.d. sequence of standard T-valued Brow- 
nian motions. In each case the corresponding initial distribution is uniform on 
[27rj72",27r(j -M)/2"[. Moreover, for n e N fixed the sequences {B"'')ieN are 
independent as j varies and the same is true of the sequences (Z"''')igN- 

Let W_ (resp. W"''-' , W^'^) be the coalescing system defined in terms of (i?^)^^^ 
(resp. (i?"'"')igN7 (■^i"'"')ieN) in the same manner that W is defined in terms of 

(■Z^OiGN- 

It is clear by construction that 



(9.16) N{t) = \W{t)\ < ^ \W'\t)\ < ^ |W:"'*(i)L t>0,ne'N. 

4=0 8=0 

Lemma 9.6. The expectation F^'^^[\W{1)\] is finite. 
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Proof. There is an obvious analogue of the duality relation Proposition 9.1 for 
systems of coalescing and annihilating one-dimensional Brownian motions. Using 
this duality and arguing as in the proof of Corollary it is easy to see that, 
letting L and U be two independent, standard, real-valued Brownian motions on 
some probability space (f2,.F, P) with Z(0) = C7(0) = 0, 



pM[|VK(l)|] 

CXD 

= lim y p[°°] {W{1) n [2TTi/M, 2TT{i + 1)/M] 7^ 0} 

i— — oo 
i—--oc 

[L{l) + 2TTi/M,U{l) + 2TT{i + l)/M]n [0,27r] 7^0 



< lim sup c' MI 



1 < mm^ {U{t) - Lit)) > -2n/M \ (L/(1) - Z(l) + c") 



for suitable constants c' and c". Noting that (f7 — L)/\/2 is a standard Brownian 
motion, the result follows from a straightforward calculation with the joint distri- 
bution of the minimum up to time 1 and value at time 1 of such a process (see, for 
example. Corollary 30 in Section 1.3 of Fre83| ). □ 



Proof of Proposition 9.4. By the Cauchy-Schwarz inequality, 

2 

\2 



(9.17) 



'N{t) \ N{t) 

1= I E^^W <Nit)J2F,itr 



i=l 



and hence, by Lemma 9.5 
(9.18) 



3 

limmf i5 7V(i) > 



On the other hand, for each n e N, \W"-'{2~^'^)\, i = 0, ... ,2" - I, are i.i.d. 
random variables which, by Brownian scaling, have the same distribution as |VF(1)|. 
By (|9J6| ), 



(9.19) 



i=0 



for 2 ^" < < < 2 ^\ An application of Lemma 9.6 and the following strong 
law of large numbers for triangular arrays completes the proof. 

Lemma 9.7. Consider a triangular array {X^.i : 1 < « < 2", n G N} of identically 
distributed, real-valued, mean zero, random variables on some probability space 
(f2,JF, P) such that the collection {X^^i : 1 < * < 2"} is independent for each 
71 e N. Then 



lim 2-" (X„4 -)-••• + X„,2.) =0, 



a.s. 
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Proof. This sort of result appears to be known in the theory of complete conver- 



gence. For example, it follows from the much more general Theorem A in | AK80| by 
taking 7V„ = 2" and ip(t) = 2* in the notation of that result (see also the Example 
following that result). For the sake of completeness, we give a short proof that was 
pointed out to us by Michael Klass. 

Let {Yn : n G N} be an independent identically distributed sequence with the 
same common distribution as the X„ i. By the strong law of large numbers, for any 
e > the probability that |Yi + • • • + F2"| > e2" infinitely often is 0. Therefore, by 
the triangle inequality, for any £ > the probability that |y2"+i + • ■ • + i^2"+i| > 
e2" infinitely often is 0; and so, by the Borel-Cantelli lemma for sequences of 
independent events, 

(9.20) J2^{\Y2^+i + ■ ■ ■ + Y2,.+i \ > 62"^} < oo 

n 

for all £ > 0. The last sum is also 

(9.21) ^P{|X„,i + ...+X„,2..| >£2"}, 

n 

and an application of the "other half" of the Borel-Cantelli lemma for possibly 
dependent events establishes that for all £ > the probability of \Xn,i + • • • + 
Xn.2^ \ > e2" infinitely often is 0, as required. □ 

10. Finitely many pure types for circular Brownian migration 

Recall that Z and Z are standard Brownian motions on the circle T and m is 
normalised Lebesgue measure. Recall also that O is the collection of open subsets 
of T that are either empty or the union of a finite number of disjoint intervals. 

Definition 10.1. Let S° denote the subset of S consisting of i' such that there 
exists a finite set {fcj, . . . , fc^} C K (depending on v) with the property that for 
m-a.e. e e T we can take i^(e) = (5^* for some i, and, moreover, we can choose a 
version of i/ such that the sets {e £ T : i/(e) — 6k*} E O for 1 < j < A^. 

Theorem 10.2. For all ^eE, Q^{Ai e S° for all t > 0} ^ 1. 
Proof. Fix /i e S and t > Q. We will first show that 
(10.1) Q^{At G S°} = 1. 



By the same argument as in Proposition 5.1 of | Eva97 |, Q'^-a.s. there is a random 
countable set of types K* such that Xt{e) G {6k ■ k G K*} for m-a.e. e G T. We 
can also require that K* has been chosen "minimally" so that m({e E E : Xt{e) = 
6k}) > for all k E K*, Q^-a.s., and this requirement specifies K* uniquely, Q'^- 
a.s. For n G N it is clear that on the event where K* has cardinality at least n 



the dissimilarity D„(At) (recall Definition 6.1) is strictly positive Q'^-a.s. It follows 



from Theorem 6.4 and Corollary |9.3| that K* is finite Q'^-a.s. 

In order to show that a representative of the equivalence class of Xt in S may 
be defined so that {e G T : Xt{e) = 6k} G O for all k E K*, it suffices by the device 
used in the proof of Theorem |6.4| to consider the case where the probability measure 
/Lt(e) is diffuse for all e G T and to show in this case that Q^-a.s. for all k G K* the 
support of the measure l{Xt{e) = 6k) m{de) (which does not depend on the choice 
of equivalence class representative) is a connected set. For this, it in turn suffices to 
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check that if ai, 61, ci, di, 02, 62, C2, d2 are arranged in anti-clockwise order around 
T, then we have 



(10.2) 



,®4 



{de)l{ei e]ai,6i[, 62 e]a2,^2[, 63 g]ci,di[, 64 e]c2,d2[} 



)Xt{e^){dh)l{k3 ^h=k2^ h} = 0, Q^-a.s. 



or, equivalently by Remark |4.4| , 

i{zW(0) e]aM 4"(0) e]«2,fo2[, 4"(0) e]ci,di[, zW(o) e]c2,d2[} 

Write, for our fixed t > 0, 



(10.3) 



(10.4) 



inf{0 <s<t:Z^ 



1 < i < j < 4, 



for the first colhsion time of and zj*^' before time t, with our standing convention 
that inf = 00. We have PW{ry = Tm ^ 00} = for (i, j) ^ {k,£). Suppose that 
we have a reahsation with the properties 

(10.5) Z^^\o) e]ai,6i[,4-'l(0) e]a2,62[,4"(0) e]ci, zfl(O) e]c2,d2[, 



(10.6) 

and 

(10.7) 



7^ 00 implies 7^ Tfe/? for (i, j) ^ {k,t). 



In order that ^^i\t) — 72'k^) holds, we must have T12 7^ 00. From the continuity 
of the paths of circular Brownian motion and (10.7), in order that (10.5) holds it 
must then be the case that 



(10.8) 



Ti3 A Tu A T23 A r24 < A T3. 



By construction, this would imply that ^^^\t) = 7i^'(^) = 72*' (^) oi' 74*'('') — 
7^1 (t) = 7;[''l(i), contradicting ( |10.6D . Thus ( |10.3| ) holds and the proof of ( |l0.lD 
is complete. 



In order to establish the claim of the theorem, it suffices by (10.1) and the 
Markov property to consider the special case of ^ G S°. Write {fcj", . . . , /c^} C K 
for the corresponding set of types k such that ra({e G T : /^(e) = 5k}) > 0. Fix 
1 < i < N. Let G C be a closed and open set such that k* G G and k* ^ G 
for j ^ i (writing k* = {hi, ft,2, ■ • ■ ) one can take G = {{h[, . . .) £ K : h[ — 
hi, . . . , /I'j = hn} for some sufficiently large n). It suffices to show for each such G 
that if we put Yt{e) := Xt(e){G) G [0, 1], then Q'^-a.s. for alH > we can choose 
a representative of Yt G L°°(T, m) such that Yt(e) G {0, 1} for m-a.e. e G T and 
{eGT:Ytie) = l}eO. 

By the remarks at the end of Section 4 of [ ]Eva97 |, we have that y is a Feller 
process with state-space the subset L°°(T, m; [0, 1]) of L°°{T, m) consisting of [0, 1]- 
valued functions (where L°°{T, m; [0, 1]) is equipped with the relative weak* topol- 



ogy). Put B := {e G T : /i(e) = 4-} G O. By the definition of X in Theorem |1 
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and Proposition 9.1, for ifj G L^(ra), 

n 

i«"(de)V®"(e)n>^t(e^) 



(10.9) 



= / m«"(de)P 



"(de)P [^-^"(6)1 jp^/'^i-- c b| 
m®"(de)V®"(e)l {{ei, . . . , e„} C V'^it)} 

n 

m«"(de)V'^"(e)niy-(t)(e. 



(recall that is defined on the probability space {T,,Q,Q)). Consequently, the 
L°°{T, to; [0, l])-valued processes Y and lys have the same finite-dimensional dis- 
tributions. Clearly, t t-^ Iv^(t) is continuous (in the weak* topology). Therefore, 
choosing our representative of Yt to be lyB(t) for alH > establishes the desired 
conclusion. □ 



11. The tree associated with coalescing circular Brownian motions 



Recall that Z and Z are standard Brownian motions on T and m is normalised 
Lebesgue measure. 



Definition 11.1. Given i,j G N, let := 

time that i and j belong to the same block. 



: M{t >0:i 
By Remark p 



j} denote the first 



the T, 



are identically 

distributed. Metrise N with the (random) metric p given by p{i,j) '■= Tij. Observe 
that p is an ultrametric; that is, the strong triangle inequality p{i, j) < p{i, k)\/p{k, j) 
holds for all i,j,k. Let (F,p) denote the completion of (N, p) . The space (F, p) is 
also ultrametric. We refer the reader to Sections 18 and 19 of | Sch84 for basic facts 
about ultrametric spaces. 



Some discussion of the space (N, p) can be found in Section 4 of [Ald93|. The 



analogue of (F, p) for another process of coale scing e xchangeable partitions of N, 
name ly Ki ngman's coalescent, is considered in | Eva98 and the counterpart of The- 
orem 11.2 below is obtained. 

For background on Hausdorff and packing dimension see [ Mat95 |. In order to 
establish some notation, we quickly recall the definitions of energy and capacity. 
Let {T,p) be a metric space. Write Afi(T) for the collection of (Borel) probability 
measures on T. A gauge is a continuous, non-increasing function / : [0, oo[— )■ [0, cxd], 
such that /(r) < oo for r > 0, /(O) = oo, and lim^^oo /(r) = 0. Given p £ Mi{T) 
and a gauge /, the energy of p, in the gauge f is the quantity 



p{dx) / p{dy) f{p(x,y)). 



The capacity of T in the gauge f is the quantity 

Cap^(T) := (inf{£/(M) : p e M^{T)})-' 
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(note by our assumptions on / that we need only consider diffuse /i € Mi{T) in the 
infimum) . 

Let Ci C [0, 1] denote the middle-^ Cantor set equipped with the usual Eu- 
clidean metric inherited f rom [0 , 1]. One of the a s sertions of the following result 
is, in the terminology of [lPP95t| (see, also, ||BP92| , |PPS96| , |Pcr96i ), that F is a.s. 
capacity-equivalent to Ci . Hence, by the results of [PPS96|, F is also a.s. capacity- 
equivalent to the zero set of (one-dimensional) Brownian motion. 



Theorem 11.2. With -probability one, the ultrametric space (F,p) is compact 
with Hausdorff and packing dimensions both equal to i . There exist random vari- 

almost surely < K* < K** < oo and for every 



ables K* , K 
gauge f 

(11.1) 



such that 



if*Cap^(Ci) < Cap^ (F) < A'**Cap^(Ci ). 

Proof. The proof is essentially a reprise of the proof of Theorem 1.1 in pLlva98 |, 
with our Proposition 9.4 and Lemma 9.5 playing the role of the statements (2.1) 



and (2.2) in |EvaS 



□ 
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